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SECTION  I 


INTRODUCTION 


The  release  of  certain  chemical  species  into  the  upper  atmosphere  results  in  luminous  clouds 
that  display  the  resonance  electronic-vibration-rotation  spectra  of  the  chemically  reacting  species. 
Such  spectra  are  seen  in  rocket  releases  of  chemicals  for  upper  atmospheric  studies,  upon 
re-entry  into  the  atmosphere  of  artificial  satellites  and  missiles,  and  as  a  result  of  energy 
deposition  in  the  atmosphere  caused  by  nuclear  weapons  effects.  Of  particular  interest  in  this 
connection  is  the  observed  spectra  of  certain  metallic  oxides.  From  band  intensity  distributions  of 
the  spectra,  and  knowledge  of  the  f-values  for  electronic  and  vibrational  transitions,  the  local 
conditions  of  the  atmosphere  can  be  determined  (Reference  1).  Such  data  are  fundamental  for  the 
analysis  of  detection  and  discrimination  problems. 

Present  theoretical  efforts,  which  are  directed  toward  a  more  complete  and  realistic  analysis 
of  the  transport  equations  governing  atmospheric  relaxation  and  the  propagation  of  artificial 
disturbances,  require  detailed  information  of  thermal  opacities  and  LWIR  absorption  in  region  of 
temperature  and  pressure  where  both  atomic  and  molecular  effects  are  important  (References  2 
and  3).  Although  various  experimental  techniques  have  been  employed  for  both  atomic  and 
molecular  systems,  theoretical  studies  have  been  largely  confined  to  an  analysis  of  the  properties 
(bound-bound,  bound-free  and  free-free)  of  atomic  systems  (References  4  and  5).  This  has  been 
due  in  large  part  to  the  unavailability  of  reliable  wavefimctions  for  diatomic  molecular  systems, 
and  particularly  for  excited  states  or  states  of  open-shell  structure.  Only  recently  (References  6-8) 
have  reUable  procedures  been  prescribed  for  such  systems  which  have  resulted  in  the  development 
of  practical  computational  programs. 

The  application  of  these  computational  methods  to  studies  of  the  electronic  structure  and 
radiation  characteristics  of  metal  oxides  has  been  reported  for  several  of  the  lighter  systems 
(References  9-11).  A  preliminary  study  of  the  uranium/oxygen  system  has  been  reported  by 
Michels  (Reference  12)  which  identified  a  large  number  of  low-lying  molecular  states  for  both  the 
UO  and  UO'*'  systems.  Of  particular  interest  was  the  discovery  of  two  structures  for  UO"*"  that 
resulted  from  two  different  spin-couplings  of  the  uranium  valence  electrons.  These  results 
suggested  strong  LWIR  radiation  in  UO+  trising  fiom  different  electronic  state  transitions. 

An  inherent  uncertainty  in  these  preliminary  calculations  was  present,  owing  to  the  rjglect 
of  relativistic  effects  that  were  much  too  difficult  to  include  in  molecular  calculations  at  that 
period  of  time.  The  7s  valence  electron  of  uranium,  and  its  corresponding  a-  bonding  molecular 
orbital,  are  highly  relativistic  in  nature  which  results  in  a  contracted  charge  density  relative  to  that 
which  would  occur  in  lighter  molecular  systems.  The  effect  of  this  contraction  on  the  relative 
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positions  of  the  low-lying  electronic  states  of  the  uranium/oxygen  system  can  now  be  calculated 
with  some  degree  of  confidence  using  newly  developed  relativistic  computer  codes. 

Because  of  inherent  difficulties  in  the  experimental  determination  of  the  spectroscopy, 

transition  probabilities  and  LWIR  radiation  for  metal  oxide  systems  and  in  light  of  the 

aforementioned  recent  progress  in  the  calculation  of  relativistic  electronic  wavefimctions, 

especially  for  diatomic  systems,  a  technical  program  for  calculating  these  properties  was 

undertaken  for  the  Defense  Nuclear  Agency  under  Contracts  DNA001-8&-C-0032, 

DNA001-85-C-0120,  DNA001-83-C-0044  and  DNA001-82-C-0015.  The  emphasis  in  this  work 

was  on  the  ions  of  uranium  and  uranium  oxide  since  these  species  have  been  determined  to  be 

important  radiators  in  the  LWIR  region.  These  studies  indicated  that,  in  addition  to  UO'*'  and 

U02''',  the  doubly  ionized  species,  U'*'  UO'*'  and  U02‘''  ■*■,  should  also  be  considered 

» 

because  of  their  role  in  charge  neutralization  processes  and  their  potential  as  early-time  radiators. 

A  summary  of  the  results  obtained  to  date  for  the  uranium/oxygen  system  indicates  that 
both  UO'*'  and  1102 are  very  strong  LWIR  radiators  in  the  10-14 /<m  region.  In  addition,  UO'*', 
which  is  a  potential  late  time  radiator  exhibits  strong  visible  bands  and  should  be  efficiently  solar 
pumped.  A  carefiil  evaluation  of  the  low-tying  excited  states  of  U02'^,  particularly  those  charge 
transfer  states  that  arise  by  promotion  of  an  electron  from  (mainly)  ligand  MO’s  to  a  central 
uranium  atom  MO,  indicates  that  such  transitions  lie  well  outside  of  the  region  for  efficient  solar 
pumping.  However,  the  density  of  electronic  states  for  U02'*'  lying  below  8  eV  is  immense,  and 
additional  studies  are  needed  for  this  system. 

Other  studies  within  DNA  chemistry  requirements  include  an  analysis  of  dielectronic 
recombination  of  e  +  U'*'  and  e  +  UO'*'  which,  are  possible  charge  neutralization  processes, 
and  an  analysis  of  the  kinetics  of  the  N  +  O2  and  O  +  N2  reactions  in  producing  vibrationally 
hot  NO  molecules  radiating  in  the  IR.  This  report  addresses  theoretical  studies  in  ail  of  the  areas 
described  above. 

The  general  composition  of  this  report  is  as  follows.  In  Section  2,  we  present  a  description  of 
the  mathematical  methods  which  were  employed  in  this  research.  Included  in  Section  2  are 
sub-sections  which  deal  with  the  construction  of  electronic  wavefimctions,  the  calculations  of 
expectation  properties,  the  evaluation  of  molecular  transition  probabilities,  and  the  calculation  of 
electronic  wavefimctions  using  both  the  ab  initio  and  density  functional  methods.  This  is  followed 
by  Section  3  which  describes  the  inclusion  of  relativistic  effects  into  the  density  fimctional  (Xa) 
method.  The  calculated  results  and  pertinent  discussitm  are  presented  in  Sectitm  4. 
Recommendations  are  presented  in  Section  5. 
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SECTION  2 


METHOD  OF  APPROACH  -  NON  RELATIVISTIC  METHODS 

2.1  QUANTUM  MECHANICAL  CALCULATIONS. 

Central  to  these  theoretical  studies  are  the  actual  quantum-mechanical  calculations  which 
must  be  carried  out  for  the  atomic  and  molecular  species.  In  particular,  the  methodology  for 
computing  the  minimum  energy  reaction  pathways  and  testing  for  vibrational  stability  must  be 
carehiUy  analyzed.  For  added  clarity,  various  aspects  of  these  calculations  are  discussed  below. 

Much  evidence  on  diatomic  and  polyatomic  systems  indicates  the  inadequacy  of  minimum 
basis  sets  for  constructing  quantitatively  correct  molecular  wavefunctions  (References  13  and  14). 
This  means  inner-shell  and  valence-shell  orbitals  of  quantiun  numbers  appropriate  to  the  atoms 
(Is,  2s,  2p,  for  C,  N,  O;  etc.).  The  main  deficiency  of  the  minimum  basis  set  is  its  inability  to 
properly  describe  polarization  and  the  change  of  orbital  shape  for  systems  which  exhibit  large 
charge  transfer  effects.  Values  of  the  screening  parameters  for  each  orbital  can  either  be  set  from 
atomic  studies  or  optimized  in  the  molecule;  the  latter  approach  is  indicated  for  studies  of  higher 
precision.  When  high  chemical  accuracy  is  required,  as  for  detailed  studies  of  the  ground  or  a 
particular  excited  state  of  a  system,  a  more  extended  basis  must  be  used.  Double-zeta  plus 
polarization  functions  or  optimized  MO’s  are  usually  required  as  a  minimum  representation  for 
reliable  calculated  results  of  chemical  accuracy. 

The  chosen  basis  sets  give  good  results  only  when  used  in  a  maximally  flexible  manner.  This 
implies  the  construction  of  perturbation  expansions  or  the  use  of  Cl  wavefunctions  with  all  kinds 
of  possible  orbital  occupancies,  so  that  the  correlation  of  electnms  into  overaU  states  can  adjust  to 
an  optimum  form  at  each  geometrical  conformation  and  for  each  state.  Except  when  well-defined 
pairings  exist,  as  for  closed  shell  and  exchange  dominated  systems,  a  single-configuration  study 
(even  of  Hartree-Fock  quality)  will  be  inadequate. 

Proper  electrom'c  states  for  systems  composed  of  light  atoms  should  possess  definite 
eigenvalues  of  the  spin  operator  as  well  as  an  appropriate  geometrical  symmetry.  The 
geometrical  symmetry  can  be  controlled  by  the  assigmnent  of  orbitals  to  each  configuration,  but 
the  spin  state  must  be  obtained  by  a  constructive  or  projective  technique.  Formulas  have  been 
developed  (Reference  15)  for  projected  construction  of  spin  states  from  orthogonal  orbitals,  and 
programs  implementing  these  formulas  have  been  in  routine  use  at  UTRC  for  several  years.  One 
of  the  least  widely  appreciated  aspects  of  the  spin-projection  problem  is  that  the  same  set  of 
occupied  spatial  orbitals  can  sometimes  be  coupled  to  give  more  than  rnie  overall  state  of  a  given 
spin  quantum  number.  It  is  necessary  to  include  in  calculations  all  such  spin  couplings,  as  the 
optimum  coupling  will  continuously  change  with  changes  in  the  molecular  ctmformation.  This  is 
especially  important  in  describing  degenerate  or  near-degenerate  excited  electronic  states. 
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In  the  sections  below,  we  describe  the  several  mathematical  approaches  that  are  applicable 
to  calculation  of  the  electronic  structure  of  molecules,  to  the  calculation  of  potential  energy 
surfaces  for  chemically  reacting  systems,  and  to  the  subsequent  calculation  of  radiation  and 
coUisional  processes.  Since  several  different  approaches  are  indicated,  we  describe  their  expected 
areas  of  applicability. 

2.1.1  Method  of  /luCu  Calculatioa. 

Zl.1.1  Bom-Oppenheimer  Separation.  For  a  system  of  n  electrons  and  N  nuclei,  and 
considering  only  electrostatic  interactions  between  the  particles,  we  have  for  the  total  Hamiltonian 


*2  r  N  N 


where 


Mel--^£v?+V*l(f„,RN) 

i*l 


and  where  m^,  m^ ,  Mj,  are  the  masses  of  the  electron,  atom  and  combined  system  mass, 
respectively.  Now  since  the  ratios  m^/ma  >  and  me/Mr  are  both  small,  (2  ‘  10~^  >  S '  1(H)  we  can 
effect  a  separation  of  the  electronic  and  nuclear  coordinates  treating  the  total  wavefimctitHi  as  a 
product  of  a  nuclear  and  an  electronic  part  We  have 


V<r„,RN)  =  ZjCk(RN)V\(rn,RN) 
k 


where  ^(Fa,  Rn)  is  an  electronic  wavefimction  parametric  in  the  nuclear  coordinates  as  given  in 
Equation  (3)  and  Xk(^N)  are  nuclear  motion  wavefunctions  which  satisfy  (neglecting  terms  of  the 
order  of  nie/ma) 
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The  cross  term  in  can  be  eb'minated  by  a  proper  change  of  variables  and  Equation 

(4)  then  reduces  to  a  3N-3  dimensional  Schrodinger  equation. 

For  most  systems,  where  the  velocity  of  motion  of  the  nuclei  is  slow  relative  to  the  electron 
velocity,  this  decoupling  of  electron  and  nuclear  motion  is  valid  and  is  referred  to  as  the  adiabatic 

approximation.  Equation  (3)  thus  defines  an  electronic  eigenstate  V'kOn*  Rn)>  parametric  in  the 

nuclear  coordinates,  and  a  corresponding  eigenvalue  EicCRn)  which  is  taken  to  represent  the 
potential  energy  curve  or  surface  corresponding  to  state  k. 

2.1.1.2  Variational  Methods.  By  an  initio  method  is  meant  one  that  starts  from  a 
zero-order  Hamiltonian  which  is  exact  except  for  relativistic  and  magnetic  effects,  and  which 
involves  the  evaluation  of  electronic  energies  and  other  relevant  quantities  for  wavefunctions 
which  are  properly  antisymmetrized  in  the  coordinates  of  all  the  electrons.  For  a  system 
containing  n  electrons  and  M  nuclei,  the  zero-order  Hamiltonian  depends  parametrically  on  the 
nuclear  positions  and  is  of  the  form 


4  N  n  jn  —  m  1 

^i-i  ri-Rj  isi<j  Ri-Rj 


M 


M 


(5) 


where  Zj  and  Ri  are  the  charge  and  position  of  nucleus  i,  q  is  the  position  of  electron  j,  and  Vj^ 

is  the  Laplacian  operator  for  electnxi  j.  All  quantities  are  in  atomic  units,  i.e.  lengths  in  bohrs, 
energies  in  hartrees  (1  hartree  »  2  Itydbergs).  The  many-electron  wavefimction  consists  of  one,  or 
a  linear  combination,  V'  ~  ^  terms  of  the  form 


*  ■AOs 

i-1 


(6) 


where  each  is  a  spatial  orbital,  A  is  the  antisynunetrizing  operator,  0,  is  the  spin-projection 
operator  for  spin  quantum  number  S,  and  9^  is  a  product  of  a  and  one-electron  spin 
functions  of  magnetic  quantum  number  M,.  No  requirement  is  imposed  as  to  the  double 
occupancy  of  the  spatial  orbital  so  linear  combinations  of  the  form  given  by  Equation  (6)  can 
describe  a  completely  general  wavefunction.  The  spatial  orbitals  ipfu  may  be  whatever  basis 

orbitals  have  been  introduced,  arbitrary  linear  combinations  thereot  or  specific  linear 
combinations  determined  pursuant  to  the  particular  calculations!  method  in  use. 
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The  spatial  orbitals  <f>fa  ,  the  spin  functions  and  the  coefficients  of  different  .  if  a 
linear  combination  of  %  is  used,  may  be  explicitly  determined  by  invoking  the  variational 

principle.  Various  specific  methods  are  described  below  for  determining  wavefiinctions.  However, 
we  should  first  observe  that  the  adequacy  of  an  ob  initio  calculation,  or  for  that  matter  any  energy 
calculation,  will  depend  crucially  upon  the  extent  to  which  the  wavefunction  can  be  qualitatively 
appropriate.  Some  of  the  considerations  surrounding  the  choice  of  a  wavefunction  are  the 
following: 

(i)  necessity  that  the  wavefunction  possess  sufficient  flexibility  to  be  able  to  describe 
dissociation  to  the  correct  atomic  and  molecular  fragments  as  various  intemuclear  separations  are 
increased; 

(ii)  maintenance  of  equivalent  quality  of  calculation  for  nuclear  geometries  differing  in  the 
nature  or  number  of  chemical  bonds; 

(iii)  ability  to  describe  degenerate  or  near-degenerate  electronic  states  when  they  are 
pertinent; 

(iv)  ability  to  describe  different  electronic  states  to  equivalent  accuracy  when  their 
interrelation  (e.g.,  crossing)  is  relevant,  in  particular,  ability  to  describe  ionic-valence  state  mixing; 

(v)  ability  to  represent  changes  in  the  coupling  of  electron  spins  as  bonds  are  broken  or 
reformed. 

The  foregoing  considerations  indicate  that  it  will  often  be  necessary  to  consider 
wavefunctions  with  more  than  a  minimum  number  of  singly-occupied  spatial  orbitals,  and  that 
there  will  be  many  potential  curves  or  surfaces  for  which  a  wavefunction  consisting  of  a  single  rpn 

cannot  suffice.  It  will  then  be  necessary  to  allow  mixing  of  with  different  degrees  of  orbital 

spatial  occupancy  so  as  to  obtain  smooth  transitions  from  the  occupancies  characteristic  of 
separated  atoms  or  molecules  (or  ions)  to  those  characteristic  of  a  compound  system  or  a 
different  fragmentation. 

Another  implication  of  the  considerations  surrounding  the  choice  of  a  wavefunction  is 
related  to  the  treatment  of  electron  spin.  Not  only  is  it  necessary  to  require  that  the  wavefunction 
be  an  eigenfunction  of  and  Sz  but  it  is  also  necessary  to  take  account  of  the  fact  that  under 
many  conditions,  there  will  be  more  than  one  spin  eigenfunction  of  given  S  and  m,.  The  different 
spin  eigenfunctions  correspond  to  different  couplings  among  the  individual  spins.  Since  reactive 
processes  involve  the  breaking  and  forming  of  electron-pair  bonds,  they  must  necessarily  be 
accompanied  by  reorganizations  of  the  spin  coupling.  A  failure  to  take  iuxount  of  this  will  lead  to 
qualitatively  inappropriate  wavefunctions. 
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In  Hartree-Fock  calculations  tp  (R)  is  restricted  to  a  single  W  which  is  assumed  to  consist 
as  nearly  as  possible  of  doubly-occupied  orbitals.  The  orbitals  are  then  selected  to  be  the 
linear  combinations  of  basis  orbitals  best  satisfying  the  variational  principle.  Writing 
-  Z  diy^y .  the  avi  are  determined  by  solving  the  matrix  Hartree-Fock  equation 


^  ~  2  Sijiavi  (each  A)  (7) 

V  V 


where  €i  is  the  orbital  energy  of  . 

The  Fock  operator  F^y  has  been  thoroughly  discussed  in  the  h'terature  (Reference  16)  and 
depends  upon  one-  and  two-electron  molecular  integrals  and  upon  the  ayi .  This  makes 
Equation  (7)  nonlinear  and  it  is  therefore  solved  iteratively.  UTRC  has  developed  programs  for 
solving  Equation  (7)  for  both  closed  and  open-shell  systems,  using  basis  sets  consisting  of  either 
Slater-type  or  Gaussian-type  atomic  orbitals.  Examples  of  their  use  are  in  the  literature 
(Reference  7). 

In  configuration  interaction  calculations,  the  overall  wavefunction  has  more  than  one  term, 
y>ft ,  and  the  Cp,  are  determined  by  invoking  the  variational  principle  to  obtain  the  secular 

equation 


2](H^-ws^)cy  *  0  (each  /t)  (g) 


where 


H^y  »  I  ^)M(R)Vy(R>/r  *  I  ’P;(R)'Py(R>/r  (9) 

Equation  (8)  is  solved  by  matrix  diagonalization  using  either  a  modified  Givens  method 
(Reference  17)  or  a  method  due  to  Shavitt  (Reference  18)  or  Raffenetti  (Reference  19). 
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The  matrix  elements  H^y  and  S^»  may  be  reduced  by  appropriate  operator  algebra  to  the 
forms 

H^y  *  (  n  ’P4?i.RN)|M(RN)P|  n  ’PviCn.Rw))  (10) 

P  i-1  i-1 

S;„,  =  (^|OsP|0M>(n'P4fi.RN)|P|  nWyi(ri,RN))  (11) 

P  i*l  i-1 

where  P  is  a  permutation  and  fp  its  pari^.  The  sum  is  over  all  permutations.  (dmjOsPldm)  is  a 

“Sanibel  coefficient”  and  the  remaining  factors  are  spatial  integrals  which  can  be  factored  into 
one-  and  two-electron  integrals.  If  the  are  orthonormal.  Equations  (10)  and  (11)  become  more 

tractable  and  the  and  S^y  may  be  evaluated  by  explicit  methods  given  in  the  literature 

(Reference  15).  Computer  programs  have  been  developed  for  carrying  out  this  procedure,  and 
they  have  been  used  for  problems  containing  up  to  106  total  electrons,  10  unpaired  electrons,  and 
several  thousand  configurations. 

The  Cl  studies  described  above  can  be  carried  out  for  any  orthonormal  set  of  for  which 
the  molecular  integrals  can  be  calculated.  Programs  developed  by  UTRC  make  specific  provision 
for  the  choice  of  the  as  Slater-type  atomic  orbitals,  as  Gaussian-type  orbitals,  as  symmetry 

molecular  orbitals,  as  Hartree-Fock  orbitals,  or  as  more  arbitrary  combinations  of  atomic 
orbitals. 

The  one-  and  two-electron  integrals  needed  for  the  above  described  method  of  calculaticm 
are  evaluated  for  STO’s  by  methods  developed  by  this  Center  (Reference  20).  For  Gaussian 
orbitals,  either  the  Camegie-Mellon  integral  package  (Reference  21)  or  the  integral  routines 
incorporated  in  the  HONDO  or  GAMESS  programs  (Reference  22)  can  be  employed. 

2.1.13  Configuration  Selection.  Using  a  double-zeta  plus  polarization  basis  set  of 
one-electron  functions,  a  typical  system  can  easily  have  of  the  order  of  10*^  configurations  in  full 
Cl  (that  resulting  from  all  possible  orbital  occupancies).  It  is  therefore  essential  to  identify  and 
use  the  configurations  describing  the  significant  part  of  the  wavefimction.  There  are  several  ways 
to  accomplish  this  objective.  First,  one  may  screen  atomic-orbital  occupancies  to  eliminate 
configurations  with  excessive  formal  charge.  Alternatively,  in  a  molecular-orbital  framework,  one 
may  eliminate  configuratimis  with  excessive  numbers  of  anti-bonding  orbitals.  A  third  possibility 
is  to  carry  out  an  initial  screening  of  configurations,  rejecting  those  whose  diagcmal  energies  and 
interaction  matrix  elements  do  not  satisfy  energy  significance  criteria. 
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Another  common  method  of  classifying  configurations  is  to  examine  the  total  number  of 
orbitals  in  the  wavefimction  that  differ  from  the  SCF  reference  wavefunction.  We  write  the  Cl 
wavefunction  as 


V'a  -  CoV'o  +  (12) 

i  j 

where  a  single  excitation  function,  <p  f  differs  from  the  SCF  reference  y>o  by  one  orbital  and  the 

double  excitation  function  (p^  by  two,  and  so  on.  The  Cl  coefficients,  Ci,  are  then  determined 

variationally  to  yield  the  lowest  possible  total  energy.  In  the  limit  of  a  complete  basis  set  (N-««o), 
and  where  all  possible  substitutions  are  included  in  tpci  >  the  variational  energy  approaches  the 

correct  nonrelativistic  Bom-Oppenheimer  result  Errors  arise  from  a  truncation  of  the  functions 
used  to  determine  the  SCF  reference  wavefunction  and  from  the  truncation  of  the  excitation 
functions  series  in  ipci  ■  The  reference  wavefunction  ipo  will  typically  be  the  same  for  a  CI  or 
Many  Body  Perturbation  Theory  (MBPT)  calculatimi;  however,  for  a  molecule  even  as  small  as 
water,  rpa  becomes  a  function  of  a  very  large  number  of  basis  functions.  Additionally,  the 

reference  state  itself  may  be  multi-dimensional  in  order  to  describe  systems  such  as  diradicals  or 
systems  with  orbital  degeneracy.  Because  of  this,  one  must  truncate  the  expansion  and  eliminate 
unimportant  configurations.  In  general  one  can  show  that 


-  {tpMP)  *  0 


(13) 


In  words,  the  matrix  elements  between  the  reference  wavefunction  and  triple  and  quadruple 
excitations  is  zero.  Thus  to  first  order,  only  single  and  double  excitations  contribute.  Because  of 
this,  many  CI  calculations  attempt  to  include  all  single  and  double  excitations  in  the  ejqpressions 
for  y;ci .  lb  go  beyond  this,  generally  more  than  one  reference  wavefunction  is  used.  Programs  to 
handle  configurations  on  all  the  above  criteria  are  available  at  UTRC. 

Other,  potentially  more  elegant  methods  of  configuration  choice  involve  formal  approaches 
based  on  natural  orbital  (Reference  23)  or  multiconfiguration  SCF  (Reference  24)  concepts,  lb 
implement  the  natural-orbital  approach,  an  initial  limited  CI  wavefunction  is  transformed  to 
natural  orbital  form,  and  the  resulting  natural  orbitals  are  used  to  form  a  new  CI.  The  desired 
result  is  a  concentration  of  the  bulk  of  the  CI  wavefunctitHi  into  a  smaller  number  of  significant 
terms.  The  multiconfiguration  SCF  approach  is  more  cumbersome,  but  in  principle  more 
effective.  It  yields  the  optimum  orbital  choice  for  a  preselected  set  of  configurations.  This 
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approach  works  well  when  a  small  number  of  dominant  configurations  can  be  readily  identified. 
The  method  is  described  briefly  below. 

2. 1.1.4  Multiconfiguration  ~  Self  Consistent  Field  Method  (MC-SCF).  The  Hartree-Fock 
self  consistent  field  method  has  been  proven  to  be  a  powerful  tool  for  the  calculation  and 
understanding  of  many  ground  state  properties  of  molecules  in  the  vicinity  of  their  equilibrium 
structure.  However,  in  most  cases  the  one  configuration  Hartree-Fock  approximation  is  not 
adequate  to  properly  describe  the  dissociation  of  molecular  bonds.  Also,  many  excited  states 
cannot  be  represented  by  a  single  configuration  wave  function.  In  order  to  calculate  properties  for 
such  states,  or  to  investigate  the  formation  of  molecular  bonds,  one  often  needs 
multiconfiguration  wave  functions  for  which  both  the  linear  coefficients,  Q  of  the  configuration 
expansion 


V  = 

i 


(14) 


as  well  as  the  set  of  orthonormal  molecular  orbitals  {^}  ,  from  which  the  cmifigurations  are 
constructed,  are  optimized  according  to  the  variational  principle.  As  is  well  known,  this 
“MC-SCF”  problem  presents  many  more  difficulties  than  the  simple  one  configuration 
Hartree-Fock  case  and  much  work  has  been  devoted  to  obtaining  convergent  solutions  during  the 
last  decade. 

The  difficulties  mainly  arise  from  the  fact  that  for  general  MC-SCF  wave  functions  the 
energy  is  not  invariant  with  respect  to  rotations  between  occupied  orbitals.  Hence,  instead  of  a 
relatively  simple  pseudo-eigenvalue  equation  in  the  one  determinant  case,  the  set  of  coupled  Fbck 
equations 


(15) 

j  i 


with  the  hermiticity  conditions 


(16) 


has  to  be  solved,  ey  are  Lagrange  multipliers  which  account  for  the  orthonormality  constraints 
imposed  on  the  orbitals.  The  Fock  operators  Fij  depend  (Hi  the  orbitals  [^]  and  the  set  of  Cl 
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coefficients  (Qj .  In  analogy  to  the  one  detenninant  case  many  attempts  have  been  made  to  solve 

these  equations  iteratively  by  keeping  the  Fock  operators  fixed  in  each  iteration  step.  Then  the 
Lagrange  multipliers  can  be  expressed  by  coupling  operators  constructed  such  that  the  Fock 
equations  are  transformed  into  pseudo-eigenvalue  equations  yielding  the  improved  orbitals.  These 
are  used  in  a  second  step  to  determine  new  Cl  coefficients  by  diagonalizing  the  Cl  matrix 

(V'i|%|V'j)  •  convergence  of  these  algorithms,  however,  has  often  been  found  to  be  poor. 

A  second  group  of  MC-SCF  methods  is  based  on  the  generalized  Brillouin  theorem.  In  these 
methods,  the  orbital  changes  are  derived  from  the  coefficients  of  a  Cl  expansion  consisting  of  the 
MC-SCF  wavefunction  and  all  one-electron  singly  excited  configurations.  A  computer  program 
(ALIS)  implementing  this  method  has  been  developed  by  Ruedenberg,  et  al.  (Reference  25).  A 
somewhat  more  elegant  program  (GAMESS)  which  also  incorporates  analytical  gradients  is  also 
available  (Reference  22).  More  recently,  various  methods  have  been  proposed  which  are  based  on 
direct  minimization  of  the  energy,  avoiding  the  Fock  operators  altogether.  Such  methods  are  now 
being  studied  in  our  laboratory  and  will  be  incorporated  into  our  existing  computer  programs  if 
they  prove  to  be  highly  efficient. 

2.1.15  Many-Body  Perturbation  Theory  (MBPT).  In  MBPT  we  again  begin  with  the  SCF 
wavefunction  (V'o)  ^  reference  and  attempt  to  account  for  the  correlation  energy.  The 
concept  of  excitation  functions  described  in  the  above  section  on  Cl  calculations  carries  over  to 
MBPT  calculations.  In  MBPT  one  can  write  the  wavefunction,  V'p  >  as 


V'p  =  e’l^o) 


(17) 


where  T  is  an  recitation  operator  defined  as 


T  =  Ts  +  Td  +  Tt  +  Tq...  (18) 

where  S,  D,  T  and  Q  refer  to  single,  double,  triple,  and  quadruple  substitutions  respectively.  One 
can  write  Tn,  in  general,  where  n  refers  to  the  number  of  excitations  as 


Tn  -  ^  2  - 

:;w 


Ijk... 

abe... 


(19) 
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where  a,b,c,...  are  excited  orbitals  and  i,j,k,...  are  orbitals  occupied  in  0o  •  The  total  energy  is  now 
given  by 


Embpt  =  (^o|3te^|0o)  (20) 

lb  evaluate  Embpt>  tf..  .  above  must  be  determined.  An  equivalent  expression  that 

makes  the  perturbation  expansion  clearer  is 


OB 


Embpt  »  2  (0o|mI(Eo-36o)'*K1M^) 

k-O 


(21) 


where  the  sum  is  over  only  so-called  linked  diagrams,  Xq  has  eigenfunctions  .  and  the 
expansion  is  of  orders  in  the  perturbation,  V  ^  (%-3^) .  The  k  >  0  term  gives  the  reference 
energy  and  for  k  >  0  correlation  corrections  are  included.  In  practice  the  MBPT  total  energy  is 
calculated  by  truncating  the  T  operator  expansion  and  projecting  onto  the  appropriate 

n-space.  This  leads  to  a  set  of  nonlinear  coupled  equations  for  the  tf  "  coefficients  which 

correspond  to  the  Cl  expansion  coefficients.  The  equations  are  solved  iteratively  and  Embpt 
evaluated.  In  practice  T4  is  an  upper  limit  that  corresponds  to  quadruple  substitutions.  The  series 
is  an  oscillatory  convergent  sum,  which  in  practice  has  proven  to  be  at  least  as  accurate  as  Eq 
with  single  and  double  excitations  included. 

The  best  possible  method  to  use  would  be  a  full  O  (all  possible  configurations)  with  a 
complete  basis  set  However,  since  the  number  of  configurations  is  proportional  to  (nj' ,  where  n 

is  the  number  of  basis  functions  and  t  is  the  le\wi  of  excitation,  it  would  be  prohibitive  to  even 
include  all  single  and  double  excitations.  This  truncation  causes  the  loss  of  size  consistency,  which 
implies  that  the  energy  calculated  for  A  and  B  as  a  molecular  system,  but  dimensionally  far  apart 
is  the  sum  of  the  energy  calculated  for  A  and  B  separately.  In  a  size  extensive  calculation  the 
energy  is  proportional  to  the  size  of  the  system.  These  properties  are  very  important  if  one  wishes 
to  compute  correct  relative  energies  on  a  potential  energy  surface,  a  necessary  criteria  for  defining 
the  reactive  pathways  of  interest  in  this  research  program.  MBPX  on  the  other  hand,  is 
guaranteed  to  have  the  correct  size-dependence  because  the  expansicMis  contain  only  the  so-called 
linked  diagrams.  In  addition,  because  of  the  computational  efficiency  of  MBPX  calculations  can 
be  performed  up  to  fourth  order  in  the  perturbation  expansion  and  include  single,  double,  triple, 
and  quadruple  excitations  in  the  calculation  of  the  correlation  energy.  Such  calculations  are 
usually  performed  with  at  least  a  split  valence  plus  polarization  basis  set  For  the  fight  element 
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compounds  which  exhibit  ionic  bonding,  the  inclusion  of  diffuse  basis  functions  and  possibly 
higher  polarization  d  or  f-functions  may  be  required  for  a  more  quantitative  treatment. 

2.1J2  Density  Functional  Approach. 

The  Xa  method  (Reference  26)  for  the  electronic  structure  of  atoms,  molecules,  clusters  and 
solids  is  a  local  potential  model  obtained  by  making  a  simple  approximation  to  the  exchange  - 
correlation  energy.  Although  this  method,  as  currently  implemented,  is  less  accurate  than 
standard  ab  initio  approaches,  it  is  very  useful  for  heavy-atom  molecular  systems.  For  such  system 
(Z  2  30) ,  relativistic  effects  become  increasingly  important  and  their  incorporation  into  density 
functional  methods  is  much  easier  than  into  ab  initio  methods.  If  we  assume  a  nonrelativistic 
Hamiltonian  with  only  electrostatic  interactions,  it  can  be  shown  that  the  total  energy  E  of  a 
system  can  be  written  exactly  (Reference  27)  (in  atomic  units)  as 


This  expression  is  exact  provided  the  Uj  are  natural  orbitals  and  nj  are  their  occupation  number 
(i.e.,  eigenfunctions  and  eigenvalues  of  the  first  order  density  matrix).  The  first  term  in  Equation 
(22)  represents  the  kinetic  and  electron-nuclear  energies.  The  second  term  is  the  nuclear  repulsion 
energy.  The  sums  (ft,  v)  are  over  all  the  nuclear  charges  in  the  system.  The  third  term  is  the 
electron-electron  repulsion  term,  which  represents  the  classical  electrostatic  energy  of  the  charge 
density  interacting  with  itself,  where 


p(l)  = 

i 


(23) 
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The  last  term  Exc  represents  the  exchange  correlation  energy  and  can  be  expressed  formally  as 


E«  =  i  I  Q(l)dTi  j (24) 

where  pxc  2)  represents  the  exchange-correlation  hole  around  an  electron  at  position  1.  In  the 
exact  expression  pxc  is  dependent  on  the  second-order  density  matrix.  In  the  Hartree-Fock 
approximation  Exc  is  the  e»:hange  energy,  ^  represents  the  Fiermi  hole  due  to  the  exclusion 

principle  and  depends  only  on  the  first-order  density  matrix.  In  the  X<,  method,  we  make  a 
simpler  assumption  about  Pxc  •  If  we  assume  that  the  exchange-correlation  hole  is  centered  on  the 

electron  and  is  spherically  symmetric,  it  can  be  shown  that  the  exchange-correlation  potential 

U«.=  f^^dra  (25) 

J  U2 

is  inversely  proportional  to  the  range  of  the  hole,  rj,  where  r*  is  defined  by 

y  r?e(l)  =  1  (26) 

Therefore,  in  the  Xq  model,  the  potential  Uxc  is  proportional  to  p^^^(r).  We  define  a  scaling 
parameter  a  such  that 


U,.(l)  -  -y(3p(l)/&r)V3  (27) 

The  expression  in  Equation  (27)  is  defined  so  that  a  ^  213  for  the  case  of  a  free  electron  gas  in 
the  Hartree-Fock  model  (Reference  28)  and  a  =>  i  for  the  potential  originally  suggested  by  Slater 
(Reference  29).  A  convenient  way  to  choose  this  parameter  for  molecular  and  solid  state 
applications  is  to  optimize  the  solutions  to  t^e  Xq  equations  in  the  atomic  limit  Schwartz 
(Reference  30)  has  done  this  for  atoms  from  z  »  1  to  z  »  41  and  found  values  between  2/3  and  1. 

In  the  “spin  polarized”  version  of  the  Xa  theory,  it  is  assumed  (as  in  the  spin-unrestricted 
Hartree-Fock  model)  that  electrons  interact  only  with  a  potential  determined  by  the  diarge 


14 


density  of  the  same  spin.  In  this  case  the  contribution  to  the  total  energy  is  summed  over  the  two 
spins,  s  =  ±  1/2. 


p,(l)U^^(l)dri 


(28) 


where  the  potential  is  spin-dependent 


(29) 


and  p,  is  the  charge  density  corresponding  to  electrons  of  spin  s.  The  spin  polarized  Xa  model 

is  useful  for  describing  atoms  and  molecules  with  open-shell  configurations  and  crystals  which  are 
ferromagnetic  or  anti-ferromagnetic. 

Once  one  has  made  the  approximation  to  the  total  energy  functional  E  in  Equation  (1), 
then  the  rest  of  the  theory  follows  from  the  application  of  the  variational  principle.  The  orbitals  Uj 
are  determined  by  demanding  that  E  be  stationary  with  respect  to  variations  in  Uj.  This  leads  to 
the  set  of  one-electron  Xa  equations 


ri2 


2ti  1 

Ui  =  €iUi 


(30) 


where  Cj  is  the  one-electron  eigenvalue  associated  with  Uj.  Since  p(r)  is  defined  in  terms  of  the 

orbitals  Uj,  Equation  (30)  must  be  solved  iteratively,  until  self-consistency  is  achieved.  Empirically, 
if  one  takes  as  an  initial  guess  that  g  is  approximately  a  sum  of  superimposed  atomic  charge 
densities,  then  the  convergence  of  this  procedure  is  fairly  rapid.  The  factor  of  2/3  multiplying  the 
potential  is  a  result  of  the  linear  dependence  of  E^  on  g .  This  also  has  a  consequence  that  the 

Xa  eigenvalues  e-,  do  not  satisfy  Koopman’s  theorem,  i.e.,  they  cannot  be  interpreted  as 
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ionization  energies.  However,  it  can  be  shown  that  the  6i  are  partial  derivatives  of  the  total 
expression  of  Equation  (22)  with  respect  to  the  occupation  number. 


dE 

drti 


(31) 


If  E  were  a  linear  function  of  nj,  then  Koopmans’  theorem  would  hold.  However,  because  of  the 
dominant  Coulomb  term,  E  is  better  approximated  by  a  quadratic  function  in  nj.  This  leads  to  the 
“transition  state”  approximation  which  allows  one  to  equate  the  difference  in  total  energy  between 
the  state  (nj,  nj)  and  (nj  -  1,  nj  +  1)  to  the  difference  in  the  one-electron  energies  €j  -  6i  calculated 
in  the  state  (nj  -  1/2,  nj  + 1/2).  The  error  in  this  approximation  is  proportional  to  third-order 
derivatives  of  E  with  respect  to  ni  and  nj,  which  are  usually  small  (Reference  31).  The  main 
advantage  of  using  the  transition  state  rather  than  directly  comparing  the  total  energy  values  is 
computational  convenience,  especially  if  the  total  energies  are  large  numbers  and  the  difference  is 
small. 

The  relationship  of  Equation  (31)  also  implies  the  existence  of  a  “Fermi  level”  for  the  ground 

state.  This  can  be  seen  by  varying  E  with  respect  to  nj  under  the  condition  that  the  sum  X  "i  is  a 

1 

constant,  i.e.. 


<3 


=  0 


(32) 


impUes  dE/dn\  =  A,  where  A  is  a  Lagrangian  multiplier.  This  implies  that  the  total  energy  is 
stationary  when  all  the  one-electron  energies  are  equal.  However,  the  occupation  numbers  are 
also  subject  to  the  restriction  0  ^  nj  ^  1.  This  leads  to  the  following  conditions  on  the  ground 
state  occupation  numbers: 

€i  <  A  •)  •  nj  =  1 
€{  >  A  *)  •  nj  *  0 
€i  =  A  •)  ■  0  s  nj  s  1 

In  other  words,  the  ground  state  eigenvalues  obey  Fermi  statistics  with  A  representing  the 
Fermi  energy.  It  should  be  noted  that,  in  contrast  to  the  Hartree-Fbck  theory,  where  all  the  nj  are 
either  0  or  1,  the  Xa  model  predicts,  in  some  cases,  fractional  occupation  numbers  at  the  Fermi 


(33) 

(34) 

(35) 
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level.  In  particular,  this  will  occur  in  a  system  (such  as  transition  metal  or  actinide  atom)  which 
has  more  than  one  open  shell. 

The  Xa  model  differs  in  other  significant  ways  from  the  Hartree-Fock  method.  In  fact,  the 
simplification  introduced  in  approximating  the  total  energy  expression  introduces  several  distinct 
advantages  over  Hartree-Fock: 

1.  The  primary  advantage  is  purely  computational.  The  one-electron  potential  in 
Equation  (30)  is  orbital-independent  and  local,  i.e.,  it  is  the  same  for  all  electrons 
(except  in  the  spin-polarized  X^,  theory)  and  is  a  multiplicative  operator.  On  the 
other  hand,  the  Hartree-Fock  potential  is  nonlocal,  or  equivalently,  there  is  a 
different  local  potential  for  each  orbital.  This  involves  a  great  deal  more 
computational  effort,  especially  for  systems  described  by  a  large  number  of  orbitals. 

It  has  been  shown  (Reference  32)  that  the  Xa  orbitals  for  the  first  and  second  row 
atoms  are  at  least  as  accurate  as  a  double-zeta  basis  set,  and  are  probably  better 
for  larger  atoms  which  involve  electrons  with  /  ^  2  . 

2  The  orbital-independent  Xq  potential  leads  to  a  better  one-electron  description 
of  electronic  excitations  of  a  system.  Both  the  unoccupied  (ni  »  0)  and  occupied 
(ni  -  1)  eigenfunctions  are  under  the  influence  of  the  same  potential  resulting  from 
the  other  N-1  electrons.  The  Hartree-Fock  virtual  orbitals  see  a  potential 
characteristic  of  the  N  occupied  orbitals,  and  therefore  are  not  as  suitable  for 
describing  the  excited  states.  Actually,  although  the  ground  state  virtual  eigenvalues 
are  usually  a  good  description  of  the  one-electron  excitations,  the  virtual  spectrum 
of  the  transition  state  potential  where  one-half  an  electron  has  been  removed  from 
the  system  gives  a  much  better  first-order  picture  of  these  levels  (Reference  33). 

3.  As  has  been  shown  by  Slater  (Reference  34),  the  Xa  model  rigorously  satisfies  both 
the  virial  and  Hellman-Feynman  theorems,  independent  of  the  value  of  the 
parameter  a  .  This  is  convenient  for  calculating  the  force  on  a  nucleus  directly  in 
terms  of  a  three-dimensional  integral,  rather  than  the  six-dimensional  integrals  in 
the  expression  for  the  total  energy  of  Equation  (22). 

Z12.1  G)mputational  Aspects  of  the  Xg  Method.  In  application  of  the  Xg  model  to  finite 
molecular  systems,  there  are  two  practical  aspects  of  the  calculations  which  must  be  considered. 
The  first  concerns  the  choice  of  the  integration  framework  for  describing  the  molecular 
wavefunctions  and  the  second  deals  with  the  choice  of  the  exchange  parameter,  a ,  in  different 
regions  of  space. 
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In  computations  with  heteronuclear  molecules,  there  are  several  free  parameters  that  must 
be  chosen:  the  ratio  of  sphere  radii  for  the  atomic  spheres  of  integration  at  a  given  intemuclear 
separation,  the  degree  of  sphere  overlap,  and  the  value  of  the  exchange  parameter  in  the  atomic 
spheres  and  the  intersphere  region. 

It  has  been  found  that  changing  the  ratio  of  the  sphere  radii  for  the  two  atoms  in  a 
heteronuclear  diatomic  molecule  introduces  changes  in  the  total  energy  that  can  be  large  on  a 
chemical  scale  (—  leV) .  A  choice  for  sphere  radii  based  on  covalent  bonding  radii  does  not 
necessarily  provide  a  good  estimate  for  these  calculations.  The  value  of  the  exchange  parameter, 
a ,  and  the  sphere  radii  and/or  sphere  overlap  is  normally  fixed  in  Xq  calculations  for  crystals 
where  the  geometry  is  fixed.  However,  to  develop  a  potential  curve,  the  molecule  description  needs 
to  change  substantially  as  the  intemuclear  separation  varies  and  the  changing  sphere  radii  include 
varying  fractions  of  the  total  molecular  charge  (Reference  35).  Studies  made  at  UTRC  have  shown 
that  at  any  given  separation  the  total  energy  calculated  from  the  Xq  model  is  a  minimum  at  the 
radii  ratio  where  the  spherically  averaged  potentials  from  the  two  atomic  centers  is  equal  at  the 
sphere  radius. 


Vi(r„)  =  V2(r,^  (36) 

This  relationship  between  the  potential  match  at  the  sphere  boimdaiy  and  the  minimum  in  the 
total  energy  appears  to  hold  exactly  for  “neutral”  atoms  and  holds  well  for  ionic  molecular 
constituents.  In  the  case  of  two  ionic  species,  the  long  range  tail  of  the  potential  must  go  like 
+  2/R  for  one  ion  and  -2/R  (in  Itydbergs)  for  the  other  ion  and  so  at  large  intemuclear 
separations,  the  tails  of  the  potential  cannot  match  well.  However,  at  reastmabie  separations,  the 
1/R  character  of  the  potential  does  not  invalidate  the  potential  match  criterion  for  radii  selection. 
This  match  for  the  atomic  potentials  is  applied  to  the  self-consistent  potentials. 

In  molecules  with  significant  charge  sharing  in  the  bonds,  the  radii  of  the  atomic  spheres  is 
frequently  increased  in  Xa  calculations  so  that  an  overlap  re^on  appears  in  the  vicinity  of  the 
bond  (Reference  36).  Studies  made  at  UTRC  show  that  the  contribution  to  the  total  molecular 
energy  from  the  exchange  integral  shows  a  minimum  at  the  optimum  sphere  radius  or  sphere 
overlap.  This  provides  a  sensitive  criterion  for  selecting  these  parameters. 

The  values  of  the  exchange  parameters  in  the  spherical  integration  region  around  each 
atomic  center  are  frequently  set  at  the  atomic  values  both  for  neutral  and  for  ionic  molecular 
constituents.  However,  for  light  atoms,  the  value  of  a  which  best  reproduces  Hartree-Fock 
results  varies  substantially  with  ionici^.  In  argon,  the  foUowing  table  compares,  for  the  neutral 
atom  and  the  positive  ion,  the  HF  energy  and  tlw  Xa  energy  calculated  for  several  values  of  a . 


18 


a 

Xa  Energy 

HF  Energy 

ArO 

.72177 

526M16 

526.8173 

Ar+V2 

.72177 

526.5857 

- 

.72213 

526.6007 

- 

Ar+^ 

.72177 

526.2447 

- 

.72213 

526.25% 

- 

.72249 

526.2745 

5262743 

The  optimum  value  of  a  changes  even  more  rapidly  in  the  fluorine  atom  going  from  0.73732  for 
to  0.72991  for  F~^.  Since  the  total  energy  depends  linearly  on  a ,  this  parameter  must  be 
chosen  carefully. 

The  intersphere  exchange  coefficient  is  chosen  to  be  a  weighted  average  of  the  atomic 
exchange  parameters  from  the  two  constituents.  At  small  intemuclear  separations,  the  optimum 
radius  for  an  atomic  sphere  frequently  places  si^ficant  amounts  of  charge  outside  that  atomic 
sphere  -  charge  that  is  still  strongly  associated  with  its  original  center  rather  than  being 
transfened  to  the  other  center  or  associated  with  the  molecular  binding  region,  lb  best  account 
for  these  cases  the  weighting  coefficients  are  chosen  to  reflect  the  origin  of  the  charge  in  the 
intersphere  (or  outersphere  region), 


«i,-Q;)  +  (a.-c® 


(37) 


where  (Q,,  -  Qf)  is  the  charge  lost  from  sphere  i  relative  to  its  atomic  value  (or  ionic  value)  Qf 
and  On  is  the  atomic  exchange  parameter  for  sphere  i.  This  value  for  ai„tenpbete  is  calculated 
dynamically  -  it  is  updated  after  each  iteration  in  the  self-consistent  calculation. 

While  for  heavy  atoms,  these  changes  in  the  exchange  parameter  would  be  small,  the  a's  for 
small  atoms  vary  rapidly  with  z  (and  with  ionici^).  The  correct  choice  of  the  exchange  parameters 
influences  not  only  the  total  energy  calculated  for  the  molecule  but  also  in  some  cases  affects  the 
distribution  of  charge  between  the  atomic  spheres  and  the  intersphere  region. 

2.U  Ab  Initio  Gaussian  \lteveAinctioD  Electronic  Structure  Codes. 

Owing  to  the  complexity  of  evaluating  multicenter  electron  repulsitm  integrals  over 
Slater-type  (exponential)  orbitals,  various  groups  have  adopted  a  computational  approach  to 
electronic  structure  calculations  based  on  gaussian  orbitals.  A  highly  devel(^)ed  series  of 
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programs,  named  GAUSSIAN  8X  (X  »  2, 6, 8),  is  available  from  Camegie-Mellon  University 
(Reference  21).  A  second  series  of  programs  has  evolved  from  the  original  version  of  Dupuis  and 
King’s  (Reference  37)  HONDO  code.  This  code  has  been  further  developed  by  separate  groups  as 
GAMESS  (Reference  22X  CADPAC  (References  38)  and  HONDO  7  (Reference  39).  Finally  a  new 
code  named  COLUMBUS  (Reference  40  and  41)  has  been  developed  by  the  Ohio 
State/ Argonne/Battelle  group.  UTRC  has  been  designated  as  a  beta  test  site  for  this  new 
development.  A  brief  description  of  the  features  of  these  codes  follows. 

2.13.1  GAUSSIAN  88/90.  GAUSSIAN  88/90  is  a  connected  system  of  programs  for 
performing  ab  initio  molecular  orbital  (MO)  calculations.  It  represents  further  development  of  the 
GAUSSIAN  70/76/80/82/86  systems  already  published.  The  contributors  to  this  program  include: 
M.  J.  Frisch,  J.  S.  Binkley,  H.  B.  Schlegel,  K.  Raghavachari,  C.  F.  Melius,  R.  L  Martin,  J.  J.  P. 
Stewart,  F.  W.  Bobrowicz,  C.  M.  Rohlfing,  L.  R.  Kahn,  D.  J.  Defrees,  R.  Seeger,  R.  A.  Whiteside, 
D.  J.  Few,  E.  M.  Fleuder,  and  J.  A.  Pople.  GAUSSIAN  8X  was  originally  implemented  on  the 
chemistry  department  DEC  VAX  11/780  computer  at  Camegie-Mellon  University.  Since  then  this 
program  has  been  installed  on  a  number  of  different  computers. 

GAUSSIAN  88/90  was  designed  with  a  transparent  input  data  stream,  making  this  program 
very  user  friendly.  All  of  the  standard  input  is  free-format  and  mnemonic.  Reasonable  defaults 
for  input  data  have  been  provided,  and  the  output  is  intended  to  be  self-explanatory.  Mechanisms 
are  available  for  the  sophisticated  user  to  override  defaults  or  interface  their  own  code  to  the 
GAUSSIAN  system.  In  this  respect,  we  have  utilized  GAUSSIAN  88/90  as  a  fundamental 
framework  for  several  applicatiems.  Options  are  being  incorporated  into  this  code  to  provide 
capabilities  beytmd  Hartree-Fock  and  various  perturbatitm  theory  options. 

The  capabilities  of  the  GAUSSIAN  88/90  system  include: 

a)  Calculation  of  one-  and  two-electron  integrals  over  s,  p,  d,  and  f  contracted  gaussian 
functions.  The  basis  functions  can  either  be  cartesian  gaussians  or  pure  angular 
momentum  functions  and  a  variety  of  basis  sets  are  stored  in  the  program  and  can  be 
requested  by  name. 

b)  Self-consistent  field  calculations  for  restricted  closed-shell  (RHF),  unrestricted 
open-shell  (UHF),  and  open-shell  restricted  (ROHF)  Hartree-Fock  wavefunctions  as 
well  as  those  types  of  multiconfigurational  wavefunctions  that  fall  within  the 
Generalized  Valence  Bond-Pftrfect  Pairing  (GVB-PP)  formalism. 

c)  Evaluation  of  various  one-electron  prO|.«Tties  of  the  Hartree-Fock  wavefunction, 
including  Mulliken  population  analysis,  multipole  moments,  and  electrostatic  fields. 

d)  Automated  geometry  optimization  to  either  minima  or  saddle  points,  and  analytical  or 
numerical  differentiation  to  produce  force  constants,  ptriarizabilities,  and  dipole 
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derivatives.  This  feattire  can  be  used  to  develop  minimum  energy  reaction  paths  along 
a  complicated  many  dimensional  potential  energy  surface. 

e)  Correlation  energy  calculations  using  M<|>ller-Plesset  perturbation  theory  carried  to 
second,  third,  or  fourth  order. 

f)  Correlation  energy  calculations  using  configuration  interaction  (Cl)  with  either  all 
double  excitations  (CID)  or  all  single  and  double  excitations  (CISD). 

g)  Correlation  energy  calculations  using  coupled  cluster  theory  with  double  substitutions 
(CCD). 

h)  Correlation  energy  calculations  using  quadratic  convergence  SCF  (QCSCF).  This  is  a 
new  highly  efficient  size-consistent  method  recently  developed  by  Pople. 

i)  Analytic  computation  of  the  nuclear  coordinate  gradient  of  the  RHF,  UHF,  ROHF, 
GVB-PP,  MP2,  CID  and  RCISD  energies. 

j)  Computation  of  force  constants  (nuclear  coordinate  second  derivatives), 
polarizabilities,  hyperpolarizabilities,  dipole  derivatives,  and  polarizability  derivatives 
either  analytically  or  numerically. 

k)  Harmonic  vibrational  analysis. 

l)  Determination  of  intensities  for  vibrational  transitions  at  the  HF,  MP2,  and  Cl  levels. 

m)  lasting  the  SCF  wavefunctions  for  stability  under  release  of  constraints. 

n)  Correlated  electron  densities  and  properties. 

o)  Minimum-energy  pathway  following  from  products  to  reactants  through  a  transition 
state-intrinsic  reaction  coordinate  finder. 

2.13.2  GAMESS.  A  wide  range  of  quantum  chemical  computations  are  possible  using 
GAMESS  (Reference  22),  a  refinement  by  M.  Schmidt  and  S.  Elbert  of  the  original  HONDO 
code  (Reference  37). 

The  capabilities  of  this  code  include: 

a)  Calculations  of  RHF/UHF/ROHF/GVB-SCF  molecular  wavefunctions. 

b)  Calculations  of  multiconfiguration  SCF  (MCSCF)  wavefunctions. 

c)  Calculations  of  Cl  wavefunctions  using  the  unitary  group  method. 

d)  Optimization  of  molecular  geometries  using  an  energy  gradient  in  terms  of  Cartesian 
or  internal  coordinates. 
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e)  Searches  for  potential  energy  surface  saddle  points. 

f)  Tracing  the  intrinsic  reaction  path  from  a  saddle  point  to  reactants  or  products. 

g)  Computation  of  normal  modes  and  vibrational  frequencies. 

h)  Calculation  of  the  following  properties: 

1.  dipole,  quadrupole,  and  octupole  moments 

Z  electrostatic  potentials 

3.  electric  field  and  electric  field  gradients 

4.  electron  density  and  spin  density 

5.  Mulliken  and  Ldwdin  population  analysis 

6.  localized  orbitals  by  the  method 

7.  virial  theorem  and  energy  components. 

GAMESS  is  a  synthesis,  with  many  major  modifications,  of  several  programs.  A  large  part  of  the 
program  is  from  HONIX).  For  pure  sp  basis  sets,  the  HONDO  symmetry  and  supermatrix 
procedure  has  been  adapted  from  GAUSSIAN  76  and  GRADSCF  integrals,  both  for  the  SCF  and 
gradient  parts.  The  GVB  section  is  a  heavily  modified  version  of  GVBONE.  A  Boys  localization 
algorithm  is  implemented  from  a  heavily  modified  version  of  Streitweiser’s  QCPE  program. 

The  Cl  module  is  based  on  Brooks  and  Schaefer’s  unitary  group  program  which  was 
modified  to  run  within  GAMESS,  using  a  Davidson  eigenvector  method  written  by  S.  T  Elbert. 
The  MCSCF  module  is  a  Newton-Raphson  procedure,  developed  at  NRCC,  based  on  the  unitary 
group  Cl  package.  The  intrinsic  reaction  coordinate  pathfinder  was  written  at  North  Dakota  State 
University. 

Z133  CADPAC.  The  Cambridge  Analytic  Derivatives  Package  (CADPAC)  (Reference  38) 
is  a  group  of  programs  which  has  been  under  development  at  Cambridge  University,  UK  It 
originated  as  a  version  of  Dupuis  and  King’s  HONDO  program.  From  its  initial  state  as  an  SCF 
gradient  package  for  closed-shell  and  UHF  wavefunctions,  this  program  has  been  extensively 
modified  with  many  of  the  old  features  being  enhanced  and  many  new  features  being  added.  The 
input  data  is  now  in  free  format  with  a  ‘keyword’  system,  to  make  the  program  easier  to  use.  T  w 
integral  routines  use  essentially  the  same  methods  as  those  in  HONDO,  being  based  cm  the  Rys 
polynomial  method,  but  have  been  extended  to  cover  f-functions.  These  routines  have  also  been 
vectorized  to  take  advantage  of  the  much  greater  availability  of  supercomputers  such  as  the  Cray 
in  recent  years.  The  SCF  programs  are  by  now  a  blend  of  techniques,  but  still  contain  a  few 
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features  from  the  initial  HONDO  program,  particularly  in  the  way  symmetry  is  handled.  There 
have  been  many  modifications  to  improve  efficiency  and  the  addition  of  level-shifting  and 
damping  techniques  and  the  implementation  of  the  DUS  method  to  aid  convergence.  A  high-spin 
open-shell  SCF  program,  and  a  completely  general  open-shell  SCF  program  have  also  been 
included. 

The  gradient  routines  have  been  considerably  altered  from  those  in  the  original  HONDO, 
the  method  having  been  changed  to  one  which  is  more  efficient  and  easier  to  extend  to 
higher-order  derivatives.  These  routines  also  work  for  f-functions.  The  original  program’s 
capabilities  for  the  calculation  of  force-constants  by  numerical  differentiation  of  gradients  and  for 
the  optimization  of  geometries  remain  essentially  intact  However  there  is  now  a  choice  of  two 
optimization  algorithms  with  the  inclusion  of  Schlegel’s  method. 

In  addition  to  the  above  modifications  a  range  of  extra  facilities  were  developed.  These 
included  a  more  powerful  method  of  calculating  one-electron  properties  and  analyzing  molecular 
charge  distributions,  and  various  ‘post-SCF  stages,  beginning  with  a  4-index  transformation. 
These  new  sections  include  M<t>ller-PIesset  perturbation  theory  for  total  energies  and  molecular 
properties,  and  coupled  Hartree-Fock  calculations  of  molecular  polarizabilities.  The  polarizability 
routines  can  calculate  dynamic  properties  at  real  or  imaginaiy  frequencies,  and  obtain  dispersion 
coefficients.  It  is  also  possible  to  use  CHF  theory  to  obtain  the  perturbations  due  to  nuclear 
displacements.  These  can  be  used  to  obtain,  analytically,  all  the  dipde  and  quadrupole  moment 
derivatives  of  a  molecule.  They  also  form  part  of  the  most  important  addition  to  the  package,  the 
section  which  calculates  analytic  second  derivatives  of  the  energy.  This  is  a  powerful  technique 
whose  speed  and  accurate  represents  a  considerable  improvement  over  numerical  differentiation. 

The  capabilities  of  the  latest  CADPAC  code  (Version  4)  include; 

a)  the  evaluation  of  one-  and  two-electrcm  integrals  over  contracted  cartesian  gaussian 
basis  functions  of  type  s,  p,  d,  or  f. 

b)  SCF  calculations  for  closed  shell,  open-shell,  UHF  and  generalized  open-shell 
techniques. 

c)  calculation  of  one-electron  properties  for  these  types  of  wavefuncdon,  including  a 
distributed  multipole  analysis. 

d)  calculation  of  the  gradients  of  the  SCF  energy. 

e)  use  of  the  gradients  for  automatic  geometry  optimiz  ,don,  and  for  the  calculadon  of 
force  constants  by  numerical  differendadon.  There  is  a  choice  of  two  opdmizadon 
algorithms. 

f)  transformadon  of  the  integrals  from  the  atomic  orbital  to  the  molecular  orbital  basis. 
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g)  M<)>ller-Plesset  perturbation  theory  calculations  to  third  order  in  the  energy  and 
second  order  in  the  one-electron  properties. 

h)  coupled  Hartree-Fock  calculations  of  polarizabilities,  including  frequency  dependence, 
and  magnetizabilities. 

i)  coupled  Hartree-Fock  calculations  of  the  perturbation  due  to  nuclear  displacements. 

j)  calculation  of  the  dipole  and  quadrupole  moment  derivatives. 

k)  calculation  of  the  second  derivatives  (force  constants)  of  the  energy  by  analytic 
methods. 

l)  analytic  calculations  of  polarizability  derivatives. 

m)  calculation  of  infrared  and  Raman  intensities,  and  the  study  of  vibrational  circular 
dichroism. 

n)  calculation  of  MP2  gradients,  dipole  moment  derivatives,  polarizabilitites  and  force 
constants  using  analytic  algorithms. 

o)  spin-projected  UHF  MP2  energies 

11.3.4  HONDO.  The  HONDO  program,  (Reference  37)  originally  developed  by  Dupuis 
and  King  at  NRCC  has  recently  been  refined  and  updated  by  M.  Dupuis  at  IBM-Kingst(Hi 
(Reference  39). 

The  following  features  are  available  in  the  present  version  (HONDO  7)  of  the  program; 

a)  Single  configuration  self-consistent-field  wavefuncticms  (closed  shell  RHF,  spin 
unrestricted  UHF,  restricted  open  shell  ROHF),  generalized  valence  bond  GVB  and 
general  multiconfiguration  self-consistent-field  MCSCF  wavefunctions,  and 
configuration  interaction  Cl  wavefimctions  can  be  calculated. 

b)  The  electron  correlation  correction  to  the  energy  of  closed  sheU  RHF  wavefunctions 
can  be  calculated  by  means  of  M4>ller-Plesset  (MP)  perturbation  theory  applied  to 
second-,  third-,  and  fourth-order  (with  or  without  the  effects  of  triple  exdtatimis). 

c)  The  effective  core  potential  approximation  can  be  used. 

d)  Optimization  of  molecular  geometries  using  the  gradient  of  the  energy  with  respect  to 
nuclear  coordinates  is  possible  with  all  but  the  Cl  and  MP  wavefunctions. 

Optimization  can  be  carried  out  in  the  cartesian  space  or  in  the  internal  coordinate 
space  with  the  possibility  of  freezing  some  cartesian  or  internal  coordinates. 
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e)  The  force  constant  matrix  in  the  cartesian  space,  and  the  vibrational  spectrum, 
including  infrared  and  Raman  intensities,  can  be  calculated  with  all  but  the  Cl  and 
MP  wavefunctions. 

f)  Calculation  of  the  dipole  moment  and  polarizability  derivatives  with  respect  to  the 
nuclear  coordinates  is  possible,  for  use  with  a  previously  calculated  force  constant 
matrix.  The  force  constant  matrix  can  be  transformed  to  the  internal  coordinate  basis. 

g)  Transition  state  structures  can  be  determined  with  all  but  the  Cl  and  MP 
wavefunctions  by  taking  advantage  of  the  energy  gradients. 

h)  The  ‘Intrinsic  Reaction  Coordinate’  (IRC)  pathway  can  be  determined  with  all  but  the 
Cl  and  MP  wavefunctions  by  taking  advantage  of  the  energy  gradients. 

i)  Molecular  energies  for  several  points  on  a  potential  energy  surface  can  be  calculated  in 
a  single  run. 

j)  Non-gradient  optimization  of  basis  function  exponents  is  possible.  The  source  code 
can  be  modified  to  cany  out  optimization  of  other  non-linear  parameters,  for  example, 
contractimi  coefficients  and  geometrical  parameters. 

k)  The  following  electronic  properties  can  be  extracted  from  the  wavefunction: 

1.  dipole  moment 

2.  quadrupole  moment 

3.  Mulliten  population,  bond  order  and  valency  analyses 

4.  spin  densi^  maps 

5.  electron  density  maps 

6.  electrostatic  potential  maps 

7.  localized  orbitals  via  Boys’  method 

8.  static  dipole  polarizability 

9.  static  first  and  second  hyperpolarizabilities. 

l)  The  potential  due  to  finite  point  charges  for  a  classical  representation  of  an 
environment,  or  a  uniform  electric  field  can  be  incorporated  into  the  one-electron 
Hamiltonian. 

2.13.5  COLUMBUS.  The  COLUMBUS  code  (Reference  40  and  41)  is  a  continuing 
development  based  on  a  joint  project  at  Ohio  State  University  (I.  Shavitt),  Uni>wrsity  of  Karlsruhe 
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(R.  Ahlrichs),  and  Argonne  National  Laboratory  (R.  Shepard).  The  unique  features  of  this  code 
are  the  incorporation  of  relativistic  or  non-relativistic  core  potentials  to  permit  analysis  of  heavy 
atom  molecular  systems  and  the  inclusion  of  a  sophisticated  Cl  package  based  on  the  unitary 
group  approach  of  Shavitt  (Reference  42).  The  basic  programs  included  in  the  COLUMBUS  code 
are  as  follows. 

1.  AO  Integrals 

This  is  R.  Pitzer’s  version  of  the  integral  package  from  HONDO  with  partial  vectorization 
of  auxiliary  integral  routines.  This  package  can  handle  up  to  g-funcdons. 

2.  Integral  Iransformation 

The  transformation  algorithm  is  written  over  all  the  orbitals  which  constitute  a  shell,  using 
contraction  coefficients  defined  over  the  primitive  basis  functions.  Limited  vectorization  is 
possible  over  the  innermost  loops  for  dyadic  operations. 

3.  SCF  Gradients 

These  routines  are  based  on  the  HONDO  version  but  include  checks  on  integral  symmetry 
to  avoid  operations  over  zero  or  near  zero  elements.  At  this  time,  only  a  first  derivative 
analysis  has  been  included. 

4.  SCF  Energy 

The  SCF  routines  are  typical  of  those  found  in  the  molecular  structure  codes.  They 
incorporate  level  shifting,  damping  and  the  incorporation  of  Pulay’s  DUS  convergence 
acceleration  procedure. 

5.  MCSCF  Analysis 

The  MCSCF  package  is  due  to  R.  Shepard  and  incorporates  extensive  vectorizatitm.  The 
output  vectors  can  be  taken  directly  from  this  package  and  transferred  to  the  Cl  program. 

6.  Multireference  Direct  Cl 

This  set  of  routines  is  based  on  the  graphical  representation  of  the  unitary  group 
approach  (GUGA)  for  constructing  a  Cl  wavefunction.  This  code  differs  frmn 
conventional  Cl  in  that  a  matrix  representation  of  the  hamilttmian  is  not  explicitly 
computed.  The  unitary  group  approach  can  be  much  more  efficient  in  most  applications 
since  ~80%  of  the  hamiltonian  matrix  elements  are  usually  zero  but  must  still  be  included. 
This  leads  to  the  storage  and  diagonalization  problems  associated  with  large,  sparse 
matrices.  However,  little  work  has  been  done  on  analysis  of  excited  states  with  the  same 
molecular  symmetry  as  the  ground  state.  An  approach  that  lo(^  promising  is  to  define  an 
approximate  vector  corresponding  to  the  excited  state  wavefunction  and  to  iteratively 
improve  upon  this  solution  using  the  direct  Cl  contraction.  The  tmcertainties  in  the 
procedure  are  applications  to  situations  where  degeneracies  or  near-degeneracies  arise. 
Further  studies  of  this  case  are  in  progress. 
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2.1.4  Spin-Proiiected  Unrestricted  Hartree-Fock  Method. 

The  unrestricted  Hartree-Fock  (UHF)  method  developed  by  Pbpie  and  Nesbet 
(References  43  and  44)  yields  the  best  single-determinant  approximation  to  the  exact  wavefunction 
for  ar;  atomic  or  molecular  system.  Such  a  wavefunction  incorporates  correlation  by  allowing 
orbitals  of  different  spin  to  adjust  to  spatially  different  forms,  thus  breaking  the  symmetry 
restrictions  of  the  conventional  (RHF)  method  (References  45  and  46).  It  is  necessary,  however,  to 
project  from  such  a  wavefunction,  a  properly  antisymmetrized  spin  and  angular  momentum  state 
in  order  to  define  eigenstates  and  eigenenergies  corresponding  to  observable  spectroscopic  states. 

Let  |a),  1^)  refer  to  the  doubly-occupied  molecular  orbitals  (MO’s),  <l>a,  ■  Let  |y),  \5),  |r) 

refer  to  the  singly  occupied  MO’s.  We  assume  that  all  MO’s  have  been  subjected  to  a 
transformation  to  orthogonal  form.  Let  be  the  coefficient  associated  with  the  permutation 

Y  *—*  d  (this  permutation  is  (-1)  times  the  overlap  of  the  permuted  spin  eigenfunction  with  the 
original  spin  eigenfunction^  and  adopt  the  convention  =  -1.  Then  the  expectation  value  of 
the  nonrelativistic  Hamiltonian 


K  =  21^0+  Zvftj) 

i  i<j 


(38) 


is  given  by 


(M)  -  22(a|U|a)  +  2(lW'>  + 

a  Y  ofi 

-  +  2 

afi  Of  Of 

yd  ^  yd 


(39) 


where  i,  j  refer  to  electron  numbers,  a,  fi..  y,  d  etc.,  refer  to  MO’s  and  the  sums  run  over  all 
a,  fi.... 
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The  equations  for  determining  the  optimum  MO’s  can  be  derived  from  the  following 
variational  form 


(36)  -  2  ^  {a  1/8)  -  +  fay  (ylo))  -X^ydivld)  =0 

afi  afi  yd  . 


(40) 


These  equations  are 

(U  +  2Jc-Kc  +  Jo-^Ko)|a )  =  ^41) 

6 


where 


(U  +  2Jc-Kc  +  Jo)|y)  +  XC^|{a|VIy)<5)  =  )  (^2) 

6  fid 


l{  <>1  V|y)<S)  =  I dr0a(r)V(r,r)0XfWO 


(wiviw)  =■  j  drw'({)v(r.r)w'(r) 


j.  =  Z0|V1« 
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Jo  -  E(«ivi«) 

d 


(43) 


K.|W)-2l(^l''|W)^)  K„|W)- 2|(«|V|W)«) 

d 


Left-multiplying  Equation  (41)  by  (  y|  and  Equation  (42)  by  ( a| ,  these  equations  can  be 
combined  to  yield 
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V|«a) 

6 


Noting  that 


XCy,{(5r|V|iy)  =  XCit(dr|V|Ty). 


Equations  (41)  and  (42)  can  be  recast  as  follows.  From  Equation  (41), 

with 


M  «  2(20,  +  1)|«  ){r(«|V|r)| 
6% 


From  Equation  (42), 


I%l^)  -  Z(2C«+  l)l^)(^r|V|<S)')  =  O.M*b) 

f  f 


fi  <5 


with 


Using  Equations  (47)  and  (48)  in  (41)  and  (42),  we  find 


(U  +  2Jc-Kc  +  Jo-^Ko-iM)la)  = 


(49) 


(U  +  2Ic-Kc  +  Jo-^Ko  +  |M^-OcM0ly)=  ) 


(50) 


Finally,  Equations  (49)  and  (50)  can  be  rearranged  to  yield  the  same  Hermitian  operator  on  the 
left-hand  side  with  the  resulting  definition  of  the  one-electron  eigenenergies  for  the  closed  and 
open-shell  eigenstates: 

fa' »(a|U  +  2Jc-Kc  + Jo-jKola  )  (51) 


€/  - (y|U  +  2Jc-Kc  +  Jo-|Ko  +  |M^|y ) 


(52) 


Equations  (51)  and  (52)  represent  the  correct  one-electron  energy  expressions  for  the 
spin-projected  eigenstate. 

This  formalism  has  been  developed  to  extend  the  conventional  HF  method  to  include 
split-shell  correlation  and  proper  spin  and  symmetry  projections.  It  is  being  incorporated  into  a 
computer  program  using  gaussian-type  orbitals  (GTO’s)  as  the  elementary  basis  functions.  A 
crucial  feature  of  this  method  is  that  dissociation  always  follows  the  lowest  energy  pathway 
thereby  permitting  a  proper  description  of  bond  formation  and  breakage.  In  contrast,  RHF  or 
MC-SCF  methods  often  exhibit  improper  dissociation  character  or  exhibit  size  inconsistencies 
owing  to  correlation  energy  changes  in  going  from  molecular  geometries  to  separated 
atom-molecule  or  atom-atom  dimensions. 

2,2  TRANSITION  PROBABILITIES. 

The  electronic  and  vibrational-rotational  wavefunctions  of  a  pair  of  states  can  be  used  to 
calculate  transition  probabilities.  If  two  molecular  states  are  separated  in  energy  by  an  amount 
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AEnm  =  hcv  (h  =  Planck’s  constant,  c  =  velocity  of  light,  v  =  frequency  in  wave  numbers),  the 
semi-classical  theory  of  radiation  (References  47  and  48)  yields  for  the  probability  of  a 
spontaneous  transition  from  an  upper  state  n  to  a  lower  state  m 


^  Snm 

3  h*c^  gn 


(53) 


Here  Anm  is  the  Einstein  coefficient  for  spontaneous  transition  from  level  n  m,  gn  is  the  total 
degeneracy  factor  for  the  upper  state 


g„  =  (2-(5^A')(2S'  +  1X2J'  +  1) 


(54) 


and  Snm  is  the  total  strength  of  a  component  line  in  a  specific  state  of  polarization  and 
propagated  in  a  fixed  direction.  A  related  quantiQr  is  the  mean  radiative  lifetime  of  state  n 
defined  by 


f  =  Z  (55) 

“  m<n 


the  summation  being  over  all  lower  levels  which  offer  aUowed  connections.  The  intensity  of  the 
emitted  radiation  is 


I 


nm 


AEmnNnAimj 


(56) 


where  Nn  is  the  number  density  in  the  upper  state  n.  This  analysis  assumes  that  all  degenerate 
states  at  the  same  level  n  are  equally  populated,  which  will  be  true  for  isotropic  excitation.  The 
total  line  strength  Snm  be  written  as  the  square  of  the  transition  moment  summed  over  all 
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degenerate  components  of  the  molecular  states  n  and  m: 


Simi  -  XMi'f 


(57) 


where  j  and  i  refer  to  all  quantum  numbers  associated  collectively  with  upper  and  lower  electronic 
states,  respectively. 

In  the  Bom-Oppenheimer  approximation,  assuming  the  separability  of  electronic  and 
nuclear  motion,  the  wavefunction  for  a  diatomic  molecule  can  be  written  as 


V'VJMA  =  V''el(r;R)^R)^JMA(^a.0)  (58) 

where  ^  electronic  wavefunction  for  state  i  at  fixed  intemuclear  separation  R,  ^R] 

is  a  vibrational  wavefunction  for  level  v,  and  V'jam  refers  to  the  rotational  state  specified 

by  electronic  angular  momentum  A ,  total  angular  momentum  J  and  magnetic  quantum  number 
M.  The  representation  is  in  a  coordinate  system  related  to  a  space-fixed  system  by  the  Eulerian 

angles  {0,x,  <l>)  ■  The  transition  moment  Mji  can  be  written,  using  the  wavefunction  given  by 
Equation  (58),  as 


Mji  *  j  +  M"}v«v''J"A  "M"*erfTvdrr  (59) 

The  subscripts  e,  v  and  r  refer  to  the  electronic,  vibrational  and  rotational  wavefunctions  and 

and  are  the  electronic  and  nuclear  electric  dipole  moments,  respectively.  Integration  over  the 
electronic  wavefunction,  in  the  Bom-Oppenheimer  approximation,  causes  the  contribution  of  the 

nuclear  moment  M°  to  vanish  for  i  j.  The  electronic  dipole  moment  can  be  written 
(References  48  and  49)  in  the  form 

M*  -  -  Xefk'  -  -  {  2efk}  •  (60) 

k  k 

where  the  primed  coordinates  refer  to  the  space  fixed  system,  the  coordinates  refer  to  a 
molectiie-fixed  system  and  is  a  group  rotation  tensor  whose  elements  are  the  direction 
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cosines  related  to  the  Eulerian  rotation  angles  {d,x,  <t>) .  Using  bracket  notation.  Equations  (59) 
and  (60)  can  be  combined  to  yield  for  the  transition  moment 

j|v'|  -XeFkliv”) .{ (61) 

k 


Mi 


MtV  J  A  Wl 

iv"J"A"M’ 


=( 


The  matrix  elements  (  J'A'M'  |g)(d,;f,0)|  )  determine  the  group  selection  rules  for  an 

allowed  transition  and  have  been  evaluated  for  many  types  of  transitions  (References  50-52). 
Summing  Equation  (61)  over  the  degenerate  magnetic  quantum  numbers  M'  and  M"  we  have 
from  Equation  (57) 


C  —  cnv’l'A'  _  #J'A’  rJN' 


(62) 


where  is  the  Hbnl  -  London  factor  (References  53  and  54)  and 


ps-  -  zi{ivi-5;Miiiv")i' 


ij 


(63) 


is  the  band  strength  for  the  transition.  G}mbining  Equations  (54X  (56)  and  (62),  we  have  for  the 
intensity  of  a  single  emitting  line  from  upper  level  n: 


Anm 


3NJ  +  1) 


(64) 


where  Nj'  is  the  number  density  in  the  upper  rotational  state  J'  and  a>n  =  (2-(5oa)  (2S'  +  1)  is 
the  electronic  degeneracy.  Ikking  an  average  value  of  for  the  whole  band.  Equation  (64) 

can  be  summed  to  yield  the  total  intensity  in  the  (v',  v")  band: 


(65) 
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where  Ny'  =  ^  Nj-  is  the  total  number  density  in  the  upper  vibrational  level  v'  and  where  we 

j' 

make  use  of  the  group  summation  property 


Xrl'!):- ~  (2J- 1)  (66) 

J" 

Comparing  Equations  (56)  and  (65)  we  have  for  the  Einstein  spontaneous  transition  coefBcient  of 
the  band  (v',  v") 


Similarly,  the  lifetime  of  an  upper  vibrational  level  v'  of  state  n  can  be  written 


(67) 


(68) 

HKD  V” 

where  the  summation  runs  over  all  v"  for  each  lower  state  m.  Equation  (67)  can  be  cast  in  the 
computational  form 


ASi..(sec'^)  »  -  [AEg;;.<a.u.)] (69) 

where  and  p^»  are  in  atomic  units.  It  is  also  often  convenient  to  relate  the  transition 

probabili^  to  the  number  of  dispersion  electrons  needed  to  explain  the  emission  strength 
classically.  This  number,  the  f-number  or  oscillator  strength  for  emission,  is  given  by 


,  ,  „  mc^h^ 

fniii,V'v"  =  . 


:ASC' 


(70) 


The  inverse  process  of  absorption  is  related  to  the  above  development  through  the  Einstein 
B  coefficient  Corresponding  to  Equation  (56X  we  have  for  a  single  line  in  absorption 


Amii 


1 


K(v)dv 


h  VmnN mBmii 


(71) 
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where  K(v)  is  the  absorption  coefGcient  of  a  beam  of  photons  of  frequency  v  and 


D  —  nnv'J’A' 
“nin  **tiiv"j''A' 


3f^c  cjjTJ”  +  1) 


(72) 


is  the  Einstein  absorption  coefficient  for  a  single  line.  Summing  over  all  lines  in  the  band  (v",  v  ), 
assuming  an  average  band  frequency,  we  obtain 


i9ax 


=  Nv 


2;r 


a^^CCDni 


Pmv 


(73) 


where  Ny”  >  ^  Nj"  is  the  total  number  density  in  the  lower  vibrational  state  v".  Corresponding 

j" 

to  Equations  (69)  and  (70)  we  can  define  an  f-number  or  oscillator  strength  for  absorption  as 


f. 


nmyv  ~ 


2mAE5;;.. 

3ft2e2a>„, 


In  computational  form.  Equation  (74)  becomes 


(74) 


Iinii,v'V 


2.AEg;.<a.u.) 
3  (Um 


(75) 


where  and  p^.>  are  in  atomic  units.  Combining  Equations  (67)  and  (70)  and  comparing 

with  Equation  (74),  we  see  that  the  absorption  and  emission  f-numbers  are  related  by 


fmnyv  =  (■^)fiiin,v'v’' 
(Om 


(76) 


Some  caution  must  be  observed  in  the  use  of  f-numbers  given  either  by  Equations  (70)  or  (74) 
since  both  band  f-numbers  and  system  f-numbers  are  defined  in  the  literature.  The  confusion 
arises  from  the  several  possible  band  averaging  schemes  that  can  be  identified. 

An  integrated  absorption  coefficient  (density  corrected)  can  be  defined  from  Equation  (73) 
as 
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(77) 


Sv",v' 


Nv''Bv".v'(l  -  exp 
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where  the  exponential  factor  corrects  for  stimulated  emission.  Equation  (77)  can  be  written  in 
terms  of  the  absorption  f-number  as 


Nv 


„  ;re‘ iNv" /,  -hcn-.v’x, 
Sv'.v - T— (1-exp— r=r^)f. 
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kT 


mn.v  V' 


(78) 


Using  hc/k  =  1.43880  cm>K,  we  obtain  a  computational  formula  for  the  integrated  absorption 
coefScient  as  Sv-yCem^^atm"*)  = 


in7/273.15wNv"w,  1.4388Ovv<v(cm-0^  , 

23795  X  10"(:j^X-f^Xl- exp - - ^)  f«n.yv 


(79) 


The  total  integrated  absorptitm  is  found  from 


Stotal  “  X  2  ( 

v"  v' 

where,  under  normal  temperature  conditions,  only  the  first  few  fundamentals  and  overtones 
contribute  to  the  summations. 

The  developments  given  above  are  rigorous  for  band  systems  where  an  average  band 
frequency  can  be  meaningfully  defined.  Further  approximations,  however,  are  often  made.  For 
example,  the  electronic  component  of  the  dipote  transition  moment  can  be  defined  as 


»ii(R)  (81) 

This  quantity  is  often  a  slowly  varying  function  of  R  and  an  average  value  can  sometimes  be 
chosen.  Equation  (63)  can  then  be  written  apprcvdmately  in  factored  form  as 

Pm"  “  <lvV'  (82) 

■j 

where  qw ,  the  square  of  the  vibrational  overlap  integral,  is  called  the  Franck-Omdmi  factor. 

%  is  evaluated  at  some  mean  value  of  the  intemuclear  separation  R.  In  addition,  it  is  smnetimes 
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possible  to  account  for  a  weak  R-dependence  in  by  a  Ikylor  series  expansion  of  this  quantity 
about  some  reference  value,  Rafi  usually  referred  to  the  (0, 0)  band.  We  have. 


%  »  +  a(R - Ra^)  +  b(R -  Re^)^  +  ... ]  (83) 

Substituting  into  Equation  (S3)  and  integrating  yields 

pSIC”  “  qvV'X  1*^^  [l  +  a(Rv'»--R^)  +  b(Rvv'-R<^)^  +  ...]  (84) 

■J 


where 


(RvV'-I^»; - — 


is  the  R-centroid  for  the  transition  and 


(85) 


(Rvv  - 


(v-|(R-R.»f|v") 

(v'|v") 


is  the  R^-centroid.  Note  that  this  last  term  differs  (to  second  order)  from  the  square  of  the 
R-centroid.  An  alternate  procedure  can  be  developed  by  evaluating  Equation  (81)  at  each 
R-centroid,  Rv  v"  •  Then 


(86) 


(87) 


Equation  (87)  assumes  that  the  vibrational  wavefrmction  product  VV’V'V"  >  behaves  like  a  delta 
function  upon  integration. 


W-VV"  =  <5{R-Rvv)(v'lv")  (88) 

The  range  of  validity  of  Equation  (87)  is  therefore  questionable,  particularly  for  band  systems 
with  bad  overlap  conditions  such  as  oxygen  Schumaim-Runge.  The  range  of  validity  of  the 
R-centroid  approximation  has  been  examined  by  Frazer  (Reference  55). 
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The  final  step  in  calculating  transition  probabilities  is  the  determination  of  9bji(R) ,  the 
electronic  dipole  transition  moment,  for  the  entire  range  of  intemuclear  separations,  R,  reached  in 
the  vibrational  levels  to  be  considered.  This  can  be  expressed  in  terms  of  a  Cl  expansion  as 

%(R)  =  )  (89) 

where  and  ci,  are  coefficients  for  vii  >  respectively. 

An  analysis  similar  to  that  yielding  Equations  (10)  and  (11)  gives 


0  n 

(v^R)|M*|V»v(R))  =  Zfp(^slOsPl^s)  {  n  V'4?k.R)|M“P|  n  Mrk.R))  (90) 

p  k>l  k-1 


The  spatial  integral  in  Equation  (90)  reduces  to  one-electron  integrals  equivalent  to  overlap 
integrals,  and  the  evaluatimi  of  Equation  (90)  can  be  carried  out  with  any  of  the  standard  ab  initio 
programs  such  as  DIATOM  or  GAMESS.  Programs  for  evaluating  %(R)  in  Equation  (89)  have 

been  developed  at  UTRC  and  examples  of  their  application  have  appeared  in  the  literature 
(Reference  8). 

For  perturbed  electronic  systems,  the  transition  dipole  moment  will  have  a  strong 
R-dependence  and  R-centroid  or  other  approximatitHis  wiU  be  invalid.  A  direct  evaluatitxi  of 
Equation  (63)  would  therefore  be  required  using  the  fully-coupled  system  of  electronic  and 
vibrational  wavefunctions  to  properly  account  for  the  source  of  the  band  perturbations. 
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SECTION  3 

DISCUSSION  OF  REIATIVISTIC  METHODS 


For  heavy  atoms  (Z  2  30) ,  and  molecular  systems  built  from  heavy  atoms,  relativistic 
effects  become  increasingly  important  and  should  be  taken  into  account  in  the  calculation  of  the 
radial  wavefunctions.  The  implementation  of  relativistic  effects  into  atomic  and  molecular 
computer  codes  is  only  fairly  recent  owing  to  the  increased  complexities  introduced  in  the 
self-consistent  field  (SCF)  procedure  and  the  greatly  increased  computer  time  required  for  such 
calculations.  Compared  with  the  non-relativistic  case,  the  Dirac-Hartree-Fock  (DHF)  method 
requires  that  two  radial  functions,  Gnij.  corresponding  to  the  large  component  and  Fnij, 
corresponding  to  the  smaU  component  must  be  calculated  for  each  of  the  two  possible  j  values. 
Thus,  the  numerical  work  of  a  DHF  relativistic  treatment  is  increased  by  nearly  a  factor  of  four 
over  the  nonrelativistic  case,  exclusive  of  increased  complexities  in  evaluation  of  the  terms  of  the 
Hamiltonian.  In  view  of  this,  methods  that  have  been  developed  to  date  for  molecular  systems 
have  involved  the  use  of  model  potentials  to  represent  relativistic  effects. 

In  the  calculation  of  the  internal  energy  of  a  molecular  system  comprised  of  n  electrons  and 
N  nuclei,  and  consideimg  only  electrostatic  interactions  between  the  particles,  we  have  for  the 
total  Hamiltonian 


where 
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(91) 


(92) 


where  me  mg ,  Mj,  are  the  masses  of  the  electron,  atom  a  and  combined  system  mass, 
respectively.  Now  since  the  ratios  me/nia  and  Mx  are  both  small,  (2  x  10^  -  5  x  10^)  we  can 
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effect  a  separation  of  the  electronic  and  nuclear  coordinates  treating  the  total  wavefunction  as  a 
product  of  a  nuclear  and  an  electronic  part.  We  have 


V<fn.  Rn)  =  X  Zk(RN)Vlt(rn,  Rn)  (93) 

k 

where  Rn)  is  an  electronic  wavefunction  parametric  in  the  nuclear  coordinates  as  given  in 
Equation  (93)  and  ;qc(Rn)  nuclear  motion  wavefunctions  which  satisfy  (neglecting  terms  of  the 
order  of  me/m^) 
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X  X  Va- +  V''(r„.RN)  Zk  = 


Ifi 
a  *  B 


(94) 


The  cross  term  in  can  be  eliminated  by  a  proper  change  of  variables  and  Equation 

(94)  then  reduces  to  a  3N-3  dimensional  Schrddinger  equation. 

For  most  systems,  where  the  velocity  of  motion  of  the  nuclei  is  slow  relative  to  the  electron 
velocity,  this  decoupling  of  electronic  and  nuclear  motion  is  valid  and  is  referred  to  as  the 

adiabatic  approximation.  Equation  (93)  thus  defines  an  electronic  eigenstate  ^(jn,  Rn)  > 

parametric  in  the  nuclear  coordinates,  and  a  corresponding  eigenvalue  E^CRn)  which  is  taken  to 
represent  the  potential  energy  curve  or  surface  corresponding  to  state  k. 

In  the  usual  ab  initio  method  for  calculating  the  electronic  properties  of  a  molecular  system, 
one  starts  from  a  zero-order  Hamiltonian  that  is  exact  except  for  relativistic  and  magnetic  effects, 
and  which  involves  the  evaluation  of  electronic  energies  and  other  relevant  quantities  for 
wavefunctions  that  are  properly  antisymmetrized  in  the  coordinates  of  all  the  electrons.  For  a 
system  containing  n  electrons  and  M  nuclei,  the  zero-order  Hamiltonian  depends  parametrically 
on  the  nuclear  positions  and  is  of  the  form 


M 


2  2  iftf] 

i-lj-1  Im  lsi<j  lsi<j 


(95) 


where  Zj  and  Rj  are  the  charge  and  position  of  nucleus  i,  rj  is  the  positimi  of  electron  j,  and 
is  the  Laplacian  operator  for  electron  j.  Ail  quantities  are  in  atomic  units,  i.e.  lengths  in  bohrs, 
energies  in  hartrees  (1  hartree  »  2  Itydbergs). 
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In  addition  to  the  electrostatic  contribution,  3Ge  t  the  complete  Hamiltonian  should  contain 
additional  terms  which  correct  for  magnetic  interactions  and  relativistic  effects.  These  correction 
terms  may  be  of  importance  in  several  applications.  These  include: 

(1)  calculation  of  the  probability  of  making  a  transition  from  one  quantum  state  to 
another  in  high-momentum  collisions  such  as  those  that  can  occur  in  hot  atom 
or  heavy  atom  chemical  dynamics  expcpments; 

(2)  determination  of  the  interaction  energy  in  heavy  nuclei  systems  such  as  Cs2  and 
UO'*',  which  exhibit  open-shell  structure  on  both  nuclei  at  infinite  intemuclear 
separations; 

(3)  calculation  of  the  intermolecular  forces  between  free  radicals,  electronically 
excited  states  of  molecules  with  open-shell  structure,  and  long  molecular 
conformations  of  possible  biological  interest 

3.1  BREIT-PAUU  HAMILTONIAN. 

The  relativistic  correction  terms  to  the  usual  electrostatic  Hamiltonian  have  been  derived 

through  order  ,  where  a  is  the  fine  structure  constant  and  are  often  referred  to  as  the 
Breit-Pauli  (Reference  56)  Hamiltonian  terms.  This  Hamiltonian  has  been  derived  by  Bethe  and 
Salpeter  (Reference  57)  for  a  two-electron  system  and  has  been  generalized  to  the  many-electron 
system  by  Hirschfelder,  et  al  (Reference  58)  and  Itoh  (Reference  59).  In  the  absence  of  external 
electric  or  magnetic  fields  we  can  represent  these  correction  terms  as  follows.  Let  ^  and 

Pj  »  jAj  denote  the  operators  for  the  spin  and  linear  moment  of  electron  j,  respectively.  Then 
the  generalized  Breit-Pauli  Hamiltom'an,  correct  to  terms  of  0(a^/M) ,  can  be  written  as: 


KbP  =  Ke  +  Kll  +  »ss  +  »LS  +  Kp  +  »D 


(96) 


where  Ke  is  given  by  Equation  (95)  and  the  correction  terms  can  be  expressed  as  follows: 
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(98) 


“d  “  ?  Tz  S  -  2  2  (101) 

L«  j  lsk<j  J 

The  first  correction  term,  36ll  .  represents  the  magnetic  orbit-orbit  coupling  terms  of  the 
electrons  arising  from  the  interaction  of  the  magnetic  fields  created  by  their  motion.  The  second 
term,  Kss .  gives  the  spin-spin  magnetic  coupling  terms  which  are  often  quite  appreciable.  For  rjk 
0,  only  the  delta-function  contribution  survives  which  represents  the  Fermi-contact  spin 
interaction.  The  third  term,  Xls  ,  is  usually  the  largest  in  magnitude  and  represents  the  spin-orbit 
interaction  between  the  spin  and  magnetic  moment  of  each  electron  and  the  spin-other  orbit 
interaction,  which  represents  the  coupling  of  the  spin  of  one  electron  with  the  magnetic  moment 
of  a  different  electron.  The  term,  36p ,  corrects  for  variation  of  the  electron  mass  with  velocity  and 

the  term,  Kd  >  represents  electron  spin  terms  identified  by  Dirac  which  appear  to  have  no 
classical  analogue. 
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Aside  from  the  spin-orbit  term,  Mls  .  usually  only  the  last  term,  Kd  .  (often  called  the 
Darwin  correction  term)  and  3£p ,  the  mass-velocity  term,  are  retained  in  the  Hamiltonian, 
yielding  the  so-called  Pauli  approximation  (Reference  57). 

The  eigenfunctions  of  the  Hamiltonian  represented  by  Equation  (%)  are  four-component 
Dirac  spinors  which  may  be  expressed  as: 


/  Piik(r)  Xkmi^yf/^)  \ 

=  iiQ„k(r)  X-Ue,(l>)) 


(102) 


where  Xkm(^><P)  are  products  of  spherical  harmonics  and  Pauli  spinors  and  Pnk(r),  Qnk(r) 
represent,  respectively,  the  large  and  small  components  of  the  radial  wave  equation.  The  exact 
solution  of  the  Breit-Pauli  Hamiltonian  has  onfy  been  given  for  one-  and  to-electron  atomic 
systems  (Reference  57)  owing  to  the  complexity  of  the  operators  for  the  general  n-electron  case. 
For  a  molecular  system,  Kolos  and  Wolniewicz  (Reference  60)  have  calculated  the  relativistic 
corrections  to  H2  using  Equation  (96)  to  0(a^.  No  heavier  molecular  systems  have  been  treated 
using  the  full  Breit-Pauli  Hamiltonian. 

APPROXIMATE  TREATMENTS. 

Although  the  Breit-Pauli  Hamiltom'an  given  in  Equation  (%)  can  formally  be  employed  in  a 
molecular  system,  both  the  multiplicity  of  terms  and  the  difficulty  of  evaluation  of  the  resultant 
molecular  integrals  has  precluded  its  general  use  to  date.  For  atomic  systems,  various 
approximate  methods  of  solution,  within  a  Hartree-fock  or  multiconfiguration  Hartree-Fbck 
framework,  have  been  proposed  for  atoms  (References  61-65).  In  most  of  these  methods,  a 
restricted  Hamiltonian  which  includes  only  the  one-electron  Dirac  terms  is  usually  employed.  The 
contributions  of  the  Breit  operators  for  spin-magnetic  interactions  and  velocity  retardation  are 
then  calculated  as  first-order  perturbations  using  the  zeroth-order  Dirac  relativistic 
wavefunctions. 

An  even  more  approximate  method  for  incorporating  the  major  relativistic  effects  has  been 
proposed  by  Cowan  and  Griffin  (Reference  66).  In  this  method,  the  mass-velocity  (Kp)  and 
Darwin  (Ko)  terms,  written  in  terms  of  the  Pauli  equation  for  one-electron  atoms,  are  simply 
added  to  the  usual  nonrelativistic  Hamiltonian  operator.  In  addition,  the  spin-orbit  terms,  v ,  are 
omitted,  thereby  reducing  the  system  of  equations  to  a  single  form  representing  the  description  of 
the  major  component  wavefunction,  Pnk(r) ,  evaluated  at  the  center-of-gravity  of  the  spin-orbit 
states.  The  rationale  for  this  approximation  lies  in  the  observation  that  detailed  atomic 
calculations  using  the  complete  DHF  method  have  indicated  that,  even  for  an  atom  as  heavy  as 
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uranium,  less  than  1  percent  of  the  total  charge  is  described  by  the  small  component  radial 
wavefunctions. 

The  resulting  equations  have  the  form: 

[- ^  +  +  vi(r)  +  HUr)  +  Hb(r)]  GUr)  = 

where  the  mass-velocity  and  Darwin  terms  take  the  form 


HWr)  - 

HWr)  - +  ^(fi-vXr)] 

and  a  »  1/137.036  is  the  fine  structure  constant  The  spin-orbit  term  is  onutted  in  Equation  (103) 
and  thus  these  equations  represent  center  of  gravity  radial  functions  averaged  over  the  two 
possible  total  angular  momentum  quantum  numbers.  Equation  (103)  represents  (apart  from  the 

neglect  of  spin-orbit  effects)  the  relativistic  corrections  to  first  order  in  .  A  more  accurate 
analysis  of  heavy  atom  energy  levels  and  spectra  is  available  through  the  use  of  the  radial 

functions,  GJy(r) ,  found  from  Equation  (103X  and  a  first-order  perturbation  calculation.  Cowan 
and  Griffin  (Reference  66)  have  illustrated  the  utility  and  accuracy  of  such  an  approach. 

Recently,  Wood  and  Boring  (Reference  67)  have  adapted  this  approximate  relativistic 
method  to  the  local  exchange  problem  and  have  implemented  the  solution  of  Equation  (103) 
within  the  context  of  the  multiple  scattering  Xq  method  (Reference  22).  The  central  field 
Hamiltonian  is  modified  to  include  mass-velocity  and  Darwin  terms,  given  by  Equations  (104) 
and  (105),  in  the  sphere  surrounding  each  atomic  center.  The  intersphere  region  in  the  multiple 
scattering  approach  (constant  potential  region)  is  treated  nonrelativistically  since  charge  in  this 
region  is  far  from  a  nucleus  and  is  screened  by  the  charge  concentrated  around  the  atomic 
centers.  The  matching  conditions  for  continuity  of  the  wavefunction  at  the  sphere  boundaries 
permits  any  necessary  charge  transfer  between  the  relativistic  intra-atomic  regions  and  the 
nonrelativistic  interatomic  constant  potential  regions.  For  an  atom,  the  Wood-Boring  treatment 
reduces  to  the  Dirac-Slater  local  exchange  method,  but  with  the  neglect  of  spin-orbit  terms. 
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The  implementation  of  Equation  (103)  into  existing  nonrelativistic  multiple  scattering 
molecular  codes  is  facilitated  by  a  change  in  the  dependent  variable,  Gj,|(r) ,  to  eliminate  the  first 
derivative  of  the  wavefunction,  illustrated  in  Equation  (105).  The  usual  Numerov  method  of 
solution  can  then  be  applied  to  the  central  field  problem;  the  only  new  requirement  being  the 
numerical  tabulation  of  the  first  and  second  derivatives  of  the  potential  at  each  grid  point  in  the 
integrations.  These  derivatives  are  computed  only  once  for  each  complete  SCF  cycle  and  thus  the 
total  required  computer  time  for  a  typical  problem  is  not  significantly  increased  as  compared  with 
a  nonrelativistic  calculation.  A  complete  self-consistent  program  incorporating  this  method  has 
been  developed  at  UTRC.  Our  code  has  been  tested  repeating  calculations  for  the  U  and  Pu 
atoms  (Reference  66),  where  we  find  excellent  agreement  with  the  more  exact,  but  cumbersome, 
Dirac-Slater  calculations.  Results  for  molecular  calculations  have  recently  been  reported  by 
Boring  and  Wood  for  UF^  and  U02*'^  (References  68  and  69).  These  calculations  were  carried  out 
to  illustrate  the  shifts  in  the  valence  levels  for  such  systems  resulting  from  relativistic  effects.  The 
total  energy  was  not  of  principal  concern. 

We  have  recently  reported  (Reference  70)  the  first  all-electron  calculation  of  the  potential 
energy  curves  for  a  molecule  (Hgj )  built  from  atoms  which  exhibit  significant  relativistic  effects. 
This  study  illustrated  that  reliable  total  energies  are  obtainable  through  a  relativistic  multiple 
scattering  density  functional  treatment,  provided  care  is  taken  to  optimize  potential  match  and 
overlap  criteria  for  such  systems.  This  study  formed  the  basis  of  the  computational  scheme  that 
we  have  employed  here  for  the  uranium/oxygen  system. 

3J  EFFECTIVE  CORE  MODELS. 

It  is  well  known  from  chemical  experience  that  the  outermost  valence  electrons  contribute 
most  to  be  determining  the  chemical  properties,  especially  spectroscopic  properties,  of  molecules. 
The  core  electrons  remain  essentially  unchanged  from  their  atomic  form  except  for  intemuclear 
separations  of  the  order  of  the  charge  radii  of  the  outer  core  region  or  less,  wherein  core 
polarization  effects  may  become  important  Since  the  computational  time  required  for  ab  initio 
calculations  of  electrom'c  structure  goes  up  at  least  quadratically  with  the  number  of  electrons  in 
the  system,  there  have  been  many  attempts  to  replace  the  more  tightly  bound  core  electrons  with 
simple  one-electron  effective  potentials  (References  71-79).  Concurrent  with  elimination  of  an 
explicit  treatment  of  the  core  electrons,  a  transformation  of  the  valence  orbital  basis  is  required  to 
insure  that  the  lowest  valence  orbital  of  each  symmetry  has  a  nodeless  radial  form,  since  it  is  well 
known  that  the  lowest  energy  eigenfunction  for  a  local  potential  must  be  nodeless  (Reference  80). 
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'^ical  of  the  several  effective  core  models  that  have  been  reported  is  that  due  to  Kahn,  et 
al.  (Reference  78)  whereby  an  effective  core  potential  is  described  in  terms  of  angularly  dependent 
projection  operators  as 


ucore  ^  uco«(r)  +  ^  £  |/m  >  [U^  -  U£‘«(r)|  <  An  I  (106) 

t  m 


where  L  is  taken  at  least  as  large  as  the  highest  angular  momentum  orbital  occupied  in  the  core. 
The  term  Uff"(r)  represents  the  effective  Coulomb  and  exchange  potential  felt  by  the  valence 
electrons.  The  second  term  essentially  accounts  for  the  repulsive  potential  between  valence  and 
core  electrons  for  each  symmetry  / .  The  only  non-local  character  exhibited  by  a  potential  of  the 
form  of  Equation  (32)  arises  from  the  /  -  dependence  which  can  be  cast  in  terms  of  one-electron 

integrals  between  the  core  and  valence  orbitals.  Explicit  two-electron  terms  connecting  core  and 
valence  orbitals  are  thus  avoided  which  greatfy  simplifies  the  calculation  of  matrix  elements  of  the 
effective  Hamiltonian.  The  potential  given  by  Equation  (106)  can  be  compared  with  the 
generalized  Phillips-Kleinman  pseudo-potential  (Reference  73). 


U^  =  2(2Jc-Kc)  +  V«>  (107) 

c 

where  Jc  and  Kc  represent  the  core  orbital  Coulomb  and  exchange  operators  and  is  a 
complicated  non-local  operator  which  guarantees  core-valence  orthogonality.  Since  the  P-K  core 
orbitals  must  simultaneously  be  eigenfunctions  of  both  the  core  and  valence  Hartree-Fbck 
Hamiltonians,  in  general,  contains  complicated  two-electron  terms  and  limits  the  usefulness 
of  Equation  (107)  over  a  full  ob  initio  treatment 

The  prescription  of  Kahn  can  be  implemented  by  analytically  fitting  a  nodeless 
pseudo-orbital,  ;ifnr>  to  a  linear  combination  of  numerical  or  analytic  Hartree-Fock  orbitals 
determined  from  a  full  self-consistent  treatment  of  the  core  electrons,  using,  for  example,  the 
multiconfiguration  Hartree-Fock  code  of  Froese-Fischer  (Reference  81).  The  comptments, 

Uf^(r) ,  of  Equation  (106)  are  then  defined  implicitly  from  the  Schrodinger  equatim 

[■?■  r  ^  (108) 
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whereby 


Uf^iT)  =  «"/  +  7  +  ^  [y-2Jv,i  +  (109) 

Equation  (108)  can  be  extended  to  relativistic  systems  in  several  ways.  Kahn,  et  al. 

(Reference  82)  suggest  an  approximate  treatment  of  adding  only  the  mass-velocity  and  Darwin 
terms  to  the  usual  electrostatic  Hamiltonian  and  to  determine  approximate  HF  orbitals  in  the 
manner  prescribed  by  Cowan  and  Griffin  (Reference  66).  Equation  (108)  is  then  used  to 
determine  an  effective  Uf*"  (r)  such  that  are  the  eigenvalues  of  the  CG  approximate 
relativistic  solution  and  Xat  curve-fitted  to  the  CG-HF  orbitals.  In  this  treatment,  the  Xnt 
represent  approximate  solutions  to  the  major  component  wavefimction,  >  determined  at  the 
center-of-gravity  of  the  spin-orbit  states. 

Lee,  et  al.  (Reference  79)  adopt  a  somewhat  more  complicated  treatment  in  which  the 
spin-orbit  operator  is  added  to  the  usual  electrostatic  Hamiltonian,  in  additicm  to  the 
mass-velocity  and  Darwin  terms  retained  in  the  Cowan-Griffin  treatment  The  large  component 
eigenfunctions  of  a  fuU  Dirac-Hartree-Fock  treatment  of  the  atom,  as  given,  for  example,  by 
Desclaux  (Reference  83)  are  then  curve-fitted  in  a  manner  similar  to  the  Kahn  treatment  but 
include  the  additional  index  for  the  particular  spin-orbit  state,  Xa^  .  Use  of  these  eigenfunctions 

in  a  molecular  system  fits  more  naturally  into  a  (J-J)  coupling  scheme  whereas  the  Xat 
determined  using  the  Kahn  method  are  more  easily  represented  using  A  -  s  coupling. 

Although  these  effective  core  models  can  often  accurately  describe  an  atomic  eigenvalue 
sequence,  including  even  high-lying  electronically  excited  states  (Reference  84^  there  are  inherent 
difficulties  in  their  application  to  molecular  environments,  where  the  maximum  angular 
momentum  component  of  the  valence  shell  orbitals  may  often  exceed  the  highest  f- value 
component  retained  in  Equation  (109).  This  is  particularly  true  for  valence  orbitals  which  exhibit 
strong  changes  from  atomic  form  through  hybridization  with  higher  angular  momentum  orbitals 
or  through  the  addition  of  more  compact  polarization  terms.  In  either  case,  since  the  relativistic 
terms  are  now  all  buried  in  a  fixed  rather  than  a  dynamic  relativistic  operator,  only  static  core 
effects  are  imposed  in  determining  the  shape  of  the  valence  molecular  orbitals.  Relativistic  effects 
between  valence  electrons  and  shielding  effects  of  the  core  by  the  valence  electrons  are  therefore 
neglected  in  these  effective  core  treatments.  In  addition,  the  models  obviously  break  down 
completely  when  the  nuclei  are  brought  together  to  dimensions  such  that  core  overlap  and 
polarization  effects  become  significant.  Unfortunately,  calculations  to  date  seem  to  indicate  that 
such  effects  begin  to  set  in  for  intemuclear  separations  of  the  order  of  equilibrium  bond  lengths. 
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SECTION  4 


DISCUSSION  OF  RESULTS 


An  analysis  of  the  electronic  structure  of  UO  and  UO'*'  using  a  relativistic  formulation  has 
been  undertaken  (References  8S  and  86).  Preliminary  calculations  were  performed  for  several 
states  of  UO  and  UO'*'  and  for  the  ground  and  excited  states  of  UO2'''  ■  An  analysis  of  the 
IRAdsible/UV  radiation  characteristics  of  selected  band  systems  in  these  uranium/oxygen  systems 
has  been  carried  out  Finally,  a  data  base  for  U^,  and  U’''^  has  been  collected  for  examining 
the  role  of  dielectronic  recombination  as  a  charge  neutralization  mechanism  in  the 
uranium/oxygen  system  (Reference  87)  and  a  preliminary  analysis  of  the  kinetic  branchings  ratios 
on  N  +  O2  reactions  has  been  performed.  A  brief  summary  of  the  results  of  these  calculations 
follows. 

4.1  UO. 

Electronic  structure  calculations  were  carried  out  for  this  system  using  a  relativistic  density 
functional  formalism.  Only  a  selected  group  of  symmetries  was  studied.  Our  calculations 
indicate  that  the  lowest  symmetry  of  UO  is  derived  from  the  (A.  S)  coupled  state.  A 

vibrational  analysis  of  the  Q  =  5  ground  state  of  UO  was  carried  out  using  a 
Hulbert-Hirschfelder  (Reference  88)  fit  to  our  calculated  potential  curves.  This  fit  yields  an 

equilibrium  intemuclear  distance  of  1.89  A  and  a  fundamental  vibrational  constant  ((»«)  of  859 

cm~^.  The  spin-orbit  interaction  was  calculated  using  U*^^  atomic  splitting  parameters.  No 
explicit  two-center  effects  are  included.  Our  calculated  spectroscopic  data  are  compared  with  the 
theoretical  work  of  Krauss  and  Stevens  (Reference  89X  the  recent  rotationally  resolved 
experimental  studies  of  Heaven  and  Nicolai  (Reference  90X  the  recent  rotational  analysis  reported 
by  Kaledin,  et  al  (Reference  91),  and  estimates  based  on  experimental  data  for  similar  systems. 
The  agreement  is  well  within  the  uncertainty  of  the  calculations  or  experimental  estimates.  The 
theoretical  studies  predict  a  ^Is  or  ground  state  for  UO  whereas  the  experimental  data  of 
Heaven  and  Nicolai  (Reference  90)  and  Kaledin  (Reference  91)  suggest  a  ^14  ground  state.  Recent 
ligand  field  model  calculations  by  Dulick  (Reference  92)  predict  a  ^Ii  ground  state  for  UO,  in 
substantial  disagreement  with  previous  theory  and  experiments.  The  crystal  field  model  neglects 
overlap  and  back  charge  transfer  effects  within  its  U'''^^  model,  both  of  which  have  been  found 
to  be  of  importance  in  more  exact  studies  carried  out  within  a  full  molecular  framework.  The 
character  of  all  of  the  low-lying  multiplets  of  UO  is  similar  however,  and  this  apparent 
discrepant  should  not  affect  our  conclusions  about  optical  or  LWIR  absorptions. 

An  analysis  of  the  LWIR  enussion  from  UO  was  carried  out  based  on  a  ground  electronic 
state.  These  calculations  should  also  be  representative  of  the  LWIR  emission  from  other 
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low-lying  electronic  states  since  they  exhibit  similar  ionicities.  Our  calculated  f-numbers, 

including  fundamentals  and  overtones,  were  reported  in  DNA-TR-82-159.  These  data 

(v'  >  v")  were  given  for  the  lowest  30  vibrational  levels.  For  UO,  our  calculated  f-number  for 

the  1-0  transidon  is  4.86  x  10~^  at  A  -  12.04/1 . 

4,2  UO+. 

Detailed  searches  of  several  symmetries  of  UO'*'  were  carried  out  to  determine  the  ground 
molecular  state  of  this  system.  Our  calculations  indicate  that  the  lowest  symmetry  of  UO'*'  is 
derived  from  the  (A,  S)  coupled  state.  A  vibrational  analysis  of  the  Q  =  9/2  ground  state  of 

UO'*'  was  carried  out  using  a  Hulbert-Hirschfelder  (Reference  88)  fit  to  our  calculated  potential 
curves.  The  spin-orbit  splittings  were  derived  from  atomic  parameters  for  the  O'** ^  ion.  This  fit 

yields  an  equilibrium  intemuclear  distance  of  1.84  A  and  a  fundamental  vibrational  constant  of 
890  cm~*.  These  data  are  compared  with  other  calculated  estimates,  since  there  are  no 
experimental  data  available.  The  agreement  between  our  work  and  that  of  Krauss  and  Stevens 
(Reference  89)  is  less  satisfactory  than  in  the  case  of  UO,  but  still  well  within  the  uncertainty  of 
the  several  calculations. 

A  perturbative  treatment  for  calculating  the  densi^  of  states  in  uranium  molecules  is 
available  through  the  use  of  ligand  field  theory.  Recent  studies  by  Dulick  (Reference  92)  predict 
that  the  ground  state  of  UO'*'  has  Q  =  3/2  in  contrast  to  the  effective  core  hamiltonian 
calculations  of  Krauss  and  Stevens  (Reference  89)  and  our  relativistic  densi^  functional 
calculations,  both  of  which  predict  that  the  ground  state  has  Q  »  9/2.  The  neglect  of  important 
molecular  effects  in  the  ligand  field  model  apparently  results  in  a  bias  toward  lower  Q  values. 

An  analysis  of  the  emission  characteristics  for  the  ground  state  of  UO***  indicates  an 
oscillator  strength  for  emission  (fro)  of  5.17  x  10~^  at  A  »  113// .  A  complete  analysis  of  our 
calculated  LWIR  emission  for  UO'*'  was  given  in  DNA-TR-82-159.  Our  calculated  LWIR 
emission  for  UO'*'  is  typical  of  that  for  a  highly  ionic  metal  oxide.  SVt  predict  strong  emission 
from  the  fundamentals  of  UO***  in  the  wavelength  region  11  —  14// .  Since  this  system  exhibits 
weak  anharmonicifrr,  we  find  the  overtones  down  in  intensity  by  several  orders  of  magnitude. 
However,  the  first  occited  state  of  UO***  (^H)  Ues  at  1200  cm-^  in  our  calculations  with  a 
predicted  electronic  oscillator  strength  of  1  x  10"*  for  the  transition.  The  electronic 

and  vibrational  manifolds  for  UO'*'  are  thus  highly  overlapped  above  the  second  vibrational  level 
of  the  ground  state. 

A  sununaiy  of  the  electronic  states  of  UO^  that  have  been  studied,  both  at  NBS  (Krauss 
and  Stevens)  and  at  UTRC,  is  shown  in  Figure  1.  Ws  find  (to  date)  nine  electronic  states  lying  in 
the  region  .4-.8//  that  are  strongly  connected  to  either  the  ground  M  manifold  or  to  the  low-lying 
manifold  of  UO'*'.  The  calculated  f-numbers  for  the  bands  1^1  -•  l^K,  1*1  -♦  3*H, 
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and  1^1^  are  given  in  Tables  1, 2, 3  and  4  respectively.  We  find  a  ctmsidmble 
shift  in  for  the  3^H  state  which  offers  a  route  for  solar  excitation  fdlowed  by  IR  radiation. 
We  also  show  in  Table  4  the  IR  absorption  onresponding  to  the  tiansition.  Owing 

to  the  vertical  nature  of  the  potential  curves,  this  behaves  similar  to  a  weak  vibrational 
transition. 

The  photoabsoiption  strengths  for  tire  and  1^I-»3^I  transitions  are 

presented  in  Tables  5-7  respectively.  We  find  that  the  1^I-»2^I  oscillator  strengths  are  large 
but  that  most  of  the  absorbed  energy  will  be  teradiated  in  the  visible  since  the  tiansition  is 
nearly  vertical.  In  contrast,  the  1^1*4  transition  (Table  6)  is  an  ingiortant  route  for 
conversion  of  solar  to  IR  photons,  since  the  oscillator  strength  is  large  and  the  R«  of  the  excited 
state  is  shifted.  This  tiansition  is  not  as  strong  as  the  1^I-43^H  transitimi,  but  is  a  contributing 
factor  in  the  total  photoabsorption  profile.  The  1^I’-»3^I  (Table  7)  tiansition  reflects  the  large 
shift  in  Re  for  the  excited  state,  but  the  tiansition  monrent  is  weak  since  this  is  an  f-4s 
excitation  on  uranium. 

In  both  the  UTRC  and  NBS  studies,  a  veiy  low-lying  excited  state  is  found  at  Te  ^ 
12(X)cm*^  Since  this  state  may  be  thermally  peculated  at  shoit  time  conditions,  the  oscillator 
strengths  coupling  this  state  to  the  known  six  low-lying  excited  states  of  UO^  have  been 
examined.  Hie  calculated  oscillator  strengths  are  given  in  Tables  8-13  respectively,  for  the 
transitions;  1^H->2*I.  1^H-42^H,  1^H-»HK,  l^H-il^H,  HH-»3<H  and  1^H-»3^L  The 
transitions  to  l^K.  l^H  and  3^1,  show  strong  oscillator  strengths  owing  to  the  shifted  Re  for 
the  excited  states.  The  overall  pattern  of  solar  pumping  from  the  state  is  thus  very  similar 

to  that  from  the  ground  1^1  state. 

Since  the  density  of  electronic  states  of  UO^  is  large  above  •>  2.0eV,  we  predict  that 
relatively  strong  solar  pun^ting,  followed  by  both  LWIR  and  visible  radiation  should  occur  for 

this  system.  This  conclusion  is  similar  tt>  that  reached  by  Krauss  and  Stevens  (Reference  89} 
based  on  their  MCSCF  analysis  of  the  UO  system.  Since  several  excited  electronic  states  of 
UO^  lying  in  the  region  of  strong  solar  flux  (.4-.8  ii),  exhibit  shifted  equilibrium  intemuclear 
separation  from  that  of  the  ground  ^I  state,  efficient  conversion  of  solar  photons  to  IR  photons 
is  found  for  this  system.  At  our  present  level  of  understanding  of  the  density  of  low-lying 
states  in  UO^,  we  predict  that  solar  pumping  LWIR  radiation  should  be  approximately  three 
times  the  earthshine  contribution. 

4.3  U02*. 

An  analysis  of  U02'*' has  been  carried  out  in  D«|,  symmetry.  McGlynn  and  Smith 
(Reference  93)  have  proposed  a  ground  state  for  this  icn  based  on  simple  molecular 
orbital  arguments.  Their  set  of  one-electron  orbital  energies  is  based  on  a  maximuin  overlap 
criterion 
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that  is  empirical  in  character.  More  modem  calculations  of  the  actinide  series  atoms  suggest  that 
a  U+5[o-2j2  or  U'*'^[0"i]2  structure  should  be  the  most  stable  configuration.  An  extensive  series 
of  calculations  of  the  possible  low-lying  symmetries  of  \J02*  now  definitely  establish  that  the 
symmetry  of  the  ground  state  is  .  but  also  predict  a  very  low-lying  ^  Au  state. 

A  series  of  calculations  of  the  low-lying  electronic  states  of  U'''^  U'*’*  and  was  first 
undertaken  to  determine  the  approximate  location  of  the  electronic  states  of  the  central  ion  which 
corresponded  to  f  -»  d  transitions.  Strong  absorption  (dipole-allowed),  corresponding  to  central 
ion  promotion,  is  predicted  to  occur  only  for  wavelengths  shorter  than  —  300  nm  for  U02'^.  A 
spin  orbit  analysis  of  these  calculated  levels  yields  the  splittings  in  J-J  coupling.  The  ^<Ps/2u 
state  is  the  predicted  ground  state  multiplet. 

In  order  to  establish  that  f  -*  d  transitions  in  U02'^  be  at  these  short  wavelengths,  a  series 
of  molecular  calculations  of  the  low-lying  electronic  states  was  undertaken.  The  lowest 
dipole-allowed  f  -•  d  transitions  on  the  central  uranium  atom  corresponds  to  the 
^<l*u(5fu)  Ag(6dg)  excitation  at  A  ~  300  nm.  This  transition  lies  beyond  the  region  of  efficient 

solar  pumping  of  UO2''’.  Hay  (Reference  94)  reports  a  similar  transition  in  his  effective  potential 
studies  of  but  at  a  somewhat  shorter  wavelength  ( 250  nm). 

A  separate  molecular  spectra  for  1102’''  arises  from  charge  transfer  .«tates  formed  by 
promotion  of  an  electron  from  (mainly)  bgand  MOs  to  a  central  uranium  atom  MO.  Careful 
convergence  studies  and  the  employroent  of  a  new  integration  scheme  give  an  equilibrium  U-O 

bond  length  in  UO2'''  of  1.73  A .  This  value,  coupled  with  a  fundamental  vj  vibrational 
frequency  of  880  cm~^  is  in  good  agreement  with  Badger's  rule  (Reference  95)  correlating  bond 
lengths  and  force  constants  for  actinide  salts  and  oxo-ions. 

Using  these  more  reUable  integration  methods,  the  low-lying  charge  transfer  states  of  U02'^ 
were  re-examined.  These  charge  transfer  states,  some  of  which  would  exhibit  veiy  strong 

absorption,  begin  at  —  47,000  cm~^  [^Hu],  with  the  first  strong  dipole  allowed  transititHis 

corresponding  to  the  excitation  frrg^  10u  -»  frrg^  l^u^ .  These  transitions  lie  at 
53,0(X)  -►  57,000  cm"^  which  correspond  to  absorption  wavelengths  in  the  VUV,  well  outside  of  the 
region  for  efficient  solar  pumping.  An  analysis  of  the  LWIR  emission  from  U02'^  has  been 
carried  out  based  on  the  ground  ^Os/2u  electronic  state.  These  calculations  should  also  be 

representative  of  the  LWIR  emission  from  other  low-lying  electronic  states  since  they  all  exhibit 
similar  ionicities.  For  U02'^,  our  calculated  f-number  for  the  1-0  transition  is  1.2  x  10*^  at 
X  a  11.43 .  This  can  be  compared  with  a  calculated  f-number  of  5.2  x  10-^  at  A  >  11.3 /i  for 
UO+.  Relatively  strong  LWIR  is  therefore  predicted  for  U02'^. 
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4.4  PHOTOELECTRON  SPECTRA  OF  UO  AND  UO2+. 

A  high  temperature  (2200-2400  K)  photoelectron  spectrum  of  UO  and  U02  has  recently 
been  reported  by  Allen,  et  al.  (Reference  96).  Although  these  data  are  of  much  lower  resolution 
than  those  reported  for  UO  by  Heaven  (Reference  97)  and  by  Kaledin  (Reference  98),  they 
represent  the  first  attempt  to  measure  the  excitation  spectra  for  U02^  in  the  gas  phase.  The 
measured  spectrum  is  shown  in  Figure  2  along  with  tentative  assignments  of  several  bands.  Allen 
finds  a  broad  band  (F')  in  the  6-8  eV  region,  corresponding  to  low-lying  excitation  in  UO"*" 
beginning  at  about  1.2  eV.  This  agrees  nicely  with  the  predicted  location  by  NBS  and  UTRC  of 
several  strong  UO***  band  systems  as  illustrated  in  Figure  1. 

The  bands  labeled  B'  and  C'  can  be  assigned  to  excitations  in  UO2  beginning  at  about 
4.3  eV  Allen,  et  al.  have  assigned  these  bands  to  02p  ionization,  which  would  represent  charge 
transfer  states  in  UO2'*'.  These  bands  should  be  very  strong  if  they  correspond  to  such  charge 
transfer  transitions.  We  have  previously  calculated  that  the  f  -•  d  transition  in  U02'^  should 
begin  at  —  4  eV  and  that  the  charge  transfer  states  should  lie  above  5.8  eV.  Allowing  for  a 
possible  spin-orbit  splitting  of  —  1  eV  in  UO2'''  states,  our  predicted  charge  transfer  states  could 
lie  as  low  as  ~  4.8  eV  thus  permitting  Allen's  and  C'  bands  assignmente  as  charge  transfer 
states.  However,  the  f  -•  d  transitions  predicted  at  4.0  eV  then  cannot  be  assigned  to  the  A' 
band  since  this  featiue  should  appear  at  —  7.6  eV  in  the  photoelectron  spectrum. 

A  more  logical  assignment  is  that  the  A',  B'  and  C'  bands  ail  belong  to  f  -►  d  transitions  in 
U02^  and  that  the  (stronger)  charge  transfer  bands  appear  as  the  unassigned  D'  system.  The 
location  of  this  band  is  in  good  agreement  with  our  calculations  which  predict  the  onset  of  charge 
transfer  states  at  —  5.8  eV  Further  studies  of  these  band  assigrunents  in  U02'^  are  in  progress 
but  these  data  of  Allen,  et  al  support  our  conclusicm  that  there  are  only  weak  absorptions  in 
U02'*’  below  0.3/i  m. 

4.5  e  +  U+,e  +  U++  OIELECTRONIC  RECOMBINATION. 

In  addition  to  our  calculations  of  the  structure  and  radiative  properties  of  metal/oxygen 
species,  an  analysis  of  the  importance  of  dielectronic  recombination  processes  was  undertaken. 
Dielectronic  recombination  (DR)  is  the  process  by  which  electron  capture  from  the  continuum  to 
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a  bound  state  is  facilitated  by  collisional  excitation  of  a  previously  bound  target  electron.  The 
process  can  be  represented  by: 


e  +  •  * 


(110) 


ju+q-i]  •  •_*  •  •  +  hv 


(111) 


where  is  a  doubly  excited  state  of  the  system  at  ionization  level  ( +  q-1)  and 

is  a  stabilized  excited  state  lying  below  the  ionization  threshold.  The  process  is  shown 
schematically  in  Figure  3.  In  the  case  of  recombination  unperturbed  by  external  fields  or  the 
plasma  density,  the  rate  coefficient  can  be  calculated  as: 


3ijjt(f) 


1.  f  1  >^ii(g)Ajk  ~-«i/kT 

L.  J  j.jj. 


(112) 


where  Uji(€)  is  the  rate  of  autoionization  from  state  |  j  >  into  the  continuum  state  |  i  > ,  Ajk  is 

the  radiative  rate  from  state  |j  >  into  the  final  bound  state  |k>  and  the  sum  is  over  all  possible 
open  channels. 

An  analysis  of  the  low-lying  electronic  states  of  U®,  and  U+2  has  been  carried  out  to 
establish  a  data  base  for  examining  the  role  of  DR  in  uranium.  The  states  examined  include  all 
one-electron  promotions  of  followed  by  one  electron  capture  final  configurations  which 
correspond  either  to  high-lying  Rydberg  states  or  bound  valence  excited  states  of  the 
system.  Ikble  14  summarizes  the  low-lying  states  which  are  important  for  DR  of  e+ and 
e+  U+  There  are  two  low-lying  configurations  for  which  have  the  J-J  coupled 
representations,  [Sfs/2^  7s]  and  [Sfs/2^  6d3/2].  These  two  configurations  are  separated  by  only 
4117  cm~i  and  both  can  be  formed  by  successively  removing  two  electrons  from  the  ground  state 
neutral  configuration,  [Sfs/2^  fids/z]-  'I^o  other  U'*'  ^  representations  are  also  possible, 
[5fs/2^  7s2]  and  [Sfsr^  7s  6d3/2],  but  these  states  lie  much  higher  in  energy  and  would  have 
negligible  population,  even  in  a  10,0(X)  K  plasma. 
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In  the  case  of  e+ U'^'DR,  at  least  two  channels  are  accessible  for  low-energy  electron 
attachment.  Since  the  Itydberg  structure  lies  s  2000  cm~*  from  the  ionization  limit,  several 
more  states  are  open  DR  channels  for  electrons  with  Te  =  2000-5000  K. 

The  low-lying  excited  DR  electronic  states  of  +  are  shown  in  Tkble  14,  relative  to  both 
the  e  +  U  +  +  [5f5/2^  7s]  and  e  +  U"*’  [5f5/2^3/2]  reactant  channels.  These  data  indicate  that  the 
first  open  attachment  channel  for  dielectronic  rea}mbination  e+U+ +[5f5/2^7s]  occurs  at 
Te  —  4000  K.  However,  the  cross-section  for  capture  into  this  channel  is  very  low  since  it  involves 
a  5Ua  7p3/2  promotion  which  has  a  low  transition  probability.  The  first  reaction  channel  for 
DR  of  e + U'*’  with  a  significant  cross-section  results  from  the  7s  -•  6d3/2  transition.  This 
channel  opens  for  electron  collision  energies,  Te>5000  K,  where  electron  capture  into  a  high  lying 
Itydberg  orbital  can  occur.  However,  the  probability  of  subsequent  autoionization  is  now 
competitive  with  radiative  stabilization  and  the  dielectronic  recombination  rate  coefficient  is  of 
the  order  of  10"^^  cmVsec  or  less.  We  conclude  that  dielectronic  recombination  of  e+ U  +  is  a 
low  probability  process  and  not  rate  competitive  with  possible  ion-neutral  charge  exchange 
reactions. 

In  contrast,  dielectronic  recombination  of  e+ U'*'  has  a  much  hi^er  piobabiliqr.  As  shown 
in  Tkble  14,  there  are  two  important  low-lying  recombination  chaimels.  A  preliminary  estimate  of 
recombination  into  the  two  lowest  lying  channels  yields  the  rates  given  in  Table  15.  The  lowest 
channel  has  a  significant  DR  rate  (10“*®  -  10"i*  cmVsec)  for  low  temperature  recombination. 

This  rate  falls  off  rapidfy  at  higher  temperature  owing  to  the  T-^^  dependence.  The  second 
channel  appears  to  peak  at  --  25000  K  with  an  indicated  DR  rate  of  5  x  10~^’  cmVsec.  Thus, 
there  may  be  some  atmospheric  conditions  where  e+ U'*'  DR  should  be  considered. 

A  corresponding  DR  process  is  possible  for  e+ UO'*'.  At  the  present  time,  only  the  optically 
connected  excited  states  of  UO'*'  have  been  studied.  The  ground  molecular  orbital  configuration 
for  UO'*'  is: 


^  I[  lo^2o^3o^4o^l;r^2;r^  1 3;rl01i5] 


(113) 


where  the  last  three  molecular  orbitals  are  nearly  pure  f-orbitals  centered  on  uranium.  Excitation 
processes  that  need  t  be  examined,  all  in  a  molecular  framework,  are:  f*  -•  d^,  f*  -*  s^ 
f^  -•  sd,  f^  -*  p^  f^  sp,  and  f^  -•  pd .  Additional  excitations  frwn  the  core  2p  level  will  be 
examined  but  it  is  expected  that  these  levels  are  too  tightly  bound  for  the  electron  energy  range  of 
interest. 
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4.6  N  +  Oi  REACTIONS. 


A  study  of  the  atmospheric  reactions  O  +  N2  and  N + O2  which  yield  vibrationally  excited 
NO  molecules  has  been  initiated.  The  sources  of  NO  emission  in  the  IR  are  not  fully  understood 
and  the  detailed  kinetic  branching  of  these  reactions  is  being  studied  (References  99  and  100).  In 
particular  the  role  of  the  N[2p]  +  02(X3  Z  g“]  reaction  in  the  production  of  vibrationally  excited 
NO  molecules  is  being  examined.  The  current  NORSE  code  assumes  that  all  the  exothermicity  of 
this  reaction  goes  into  vibrational  excitation  of  the  NO  product  molecule,  with  all  levels  being 
equally  populated  up  to  v  =  26.  The  examination,  in  detail,  of  the  kinetic  routes  of  the  N + O2 
reaction,  including  those  which  lead  to  electronic  as  well  as  vibrationally  excited  products  is 
described  below.  A  preliminary  ab  initio  study  of  the  reaction  has  been  carried  out  which 
suggests  that  only  a  fraction  of  the  collisions  result  in  vibrational  excitation  in  the  product 
channels. 

The  molecular  correlation  diagram  for  the  N + O2  reaction  is  given  in  Ikble  16  for  Cj,  C2v 
and  Co.v  symmetry.  This  corresponds,  respectively,  to  oblique,  perpendicular  and  linear 
collisions  of  the  N  atom  with  the  O2  target  molecule.  The  adiabatic  molecular  correlations  for 
C«v  and  Q  symmetry  are  shown  in  Figures  4  and  S,  respectively. .  These  diagrams  represent  a 
more  complete  version  of  those  given  by  Donovan  and  Husain  (Reference  iOl)  and  by  Schofield 
(Reference  102). 

Since  the  kinetic  behavior  is  governed  by  an  average  over  all  collisional  trajectories,  we 
illustrate  the  adiabatic  correlations  in  C,  symmetry  for  the  doublet  surfaces  in  Figure  6,  and  for 
the  quartet  surfaces  in  Figure  7.  In  particular,  we  illustrate  the  adiabatic  connectitms  from 
Nl^P]  +  02[X3  Z  g~]  in  bold  lines  on  these  figures.  The  energetics  of  the  adiabatic  reaction 
surfaces  arising  from  N{2p]  +  O2  and  N^^D] + O2  are  shown  in  Ihbles  17  and  18,  respectively.  For 
N[2P]  +  O2,  only  the  third  reaction  [yielding  NOpP  IT  ]  +  0*(^D)]  is  exothermic  enough  to  yield 
NO  with  a  significant  degree  of  vibrational  excitation,  provided  that  this  reaction  proceeds  along 
an  adiabatic  pathway.  However,  there  is  spectroscopic  evidence  that  curve-crossing  and 
non-adiabatic  behavior  may  occur  for  this  system.  The  ground  ^Ai  state  of  NO2  adiabatically 
correlates  to  the  reactants:  N^S]  +  02[X3  Z  g~},  whereas  the  first  excited  state  correlates  to 
excited  state  reactants,  N[2D]  +  02[X^  Z  g-].  However,  the  ^Bi  state  lies  lower  as  we  connect  to 
the  NO[X^  n  ]  +  0[^P]  dissociation  limit.  In  addition,  Gilmore  (Reference  103)  has  pointed  out 
that  NO  quenching  of  both  0*[^S]  and  0*[iD]  proceeds  veiy  rapidly,  further  indicating 
curve-crossings  in  tlie  doublet  surfaces. 

In  order  to  clarify  the  branching  kinetics  of  the  N*[2p]  +  02(X3Zg“]  reaction,  we  have 
initiated  an  extensive  series  of  calculations  of  the  potential  energy  reaction  surfaces  for  NO2  in  the 
following  symmetries:  ^A',  ^A",  ^A',  *A".  Preliminary  Cl  calculations,  using  the  GAMESS  code 
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(Reference  104),  were  carried  out  in  C,  symmetry  assuming  a  frozen  core  of 
(la'2X2a'2X3a'2X4a'2X5a'2X6a'2),  which  is  the  Hartree-Fock  representation  of  the  inner  Is  and  2s 
electrons.  Unfortunately,  the  characteristic  charge  density  corresponding  to  this  core 
representation  differs  significantly  for  different  spin  and  spatial  couplings  of  N + 02.  The 
characteristic  charge  of  the  (6a '^)  pair,  for  example,  differs  in  calculations  in  ^A'  symmetry  as 
compared  with  ^A"  or  symmetry.  Thus  only  a  frozen  inner  shell  (Is  electrons)  appears  to  be 
possible  for  this  system.  This  smaller  frozen  core  representation  results  in  a  great  increase  in  the 
Cl  size,  from  3048  to  28S03  for  ^A'  symmetry  and  from  3000  to  2812S  for  ^A"  symmetry.  The 
calculations  must  therefore  be  carried  out  on  a  CRAY  computer  and  the  necessary  code 
modifications  for  running  this  problem  on  the  DNA  access  machines  at  LANL  have  been  made. 
The  calculations  are  carried  out  in  Q  symmetry  since  we  need  to  simultaneously  examine  the 
N+O2  reaction  path  on  nineteen  surfaces  in  ^A'  and  ^A"  symmetry  and  on  twelve  surfaces  in 
*A'  and  *A"  symmetry. 

The  results  of  our  initial  studies  of  the  long  range  N + O2  behavior  are  suirunarized  in  lable 
19.  We  find  that  two  of  the  states  correlating  to  N{2p]  +  02PC^  2  g“],  ^A'  X  and  ^A"  K,  exhibit 
long  range  repulsion  and  thus  should  correlate  adiabatically  to  NO[a^  11  ]  +  0[^P]  and 
NO[X2  n  ]  +  0[^S]  respectively.  This  leaves  <miy  the  ^A"  Vm  surface  correlating  to 
NO(X2  n  ]  +  0('D].  For  the  quartet  surfaces  which  adiabatically  correlate  to  NO[a^  11  ]  +  0[^P], 
no  vibrationally  excited  ground  state  NOfX^Il  ]  is  predicted.  Detailed  calculations  of  these 
quartet  surfaces  are  in  progress. 

The  statistical  branching  in  the  N+O2  reaction  can  be  partially  analyzed  in  terms  of  these 
calculated  potential  energy  surfaces.  In  Ihble  17,  which  lists  the  N*[2P]  +  O2  reaction  surfaces, 
we  see  that  only  the  third  reaction,  which  yields  NOpC^  11  ]  +  0[^D],  is  exothermic  enough  to 
yield  NO  with  any  significant  degree  of  vibrational  excitation.  The  reaction  surfaces  for  N*[^D] 

+  O2  are  listed  in  Ikble  18.  Three  surfaces  represent  closed  (endothermic)  channels  for  low 
energy  collisions  and  are  non-reactive  for  air  chemistry.  Six  surfaces  produce  ground  state  NO 
with  nearly  4  eV  internal  (rovibrational)  energy  and  one  surface  produces  NO  with  less  than  2  eV 
internal  energy.  Thus,  of  the  16  potential  surfaces  governing  the  Nf^D]  +  O2  and  N[2P]  +  O2 
reactions,  we  find  the  following; 

1.  3  are  non-reactive  (endothermic) 

2.  6  lead  to  electronically  excited  products  with  the  result  that  the  NO  that  is  produced 
is  not  rovibrational  excited. 

3.  7  lead  to  the  production  of  ground  state  NO  with  significant  (3-4  eV)  rovibrational 
excitation. 


Further,  60%  of  all  collisions  involving  Nf^D]  atoms  should  produce  NO  with  significant 
internal  excitation  but  only  1/6  of  all  collisions  involving  N[^P]  atoms  appear  to  lead  to 
rovibrationally  excited  NO  molecules.  Non-adiabatic  behavior  at  closer  intemuclear  separations 
could  modify  these  relative  branching  ratios  but  will  not  change  the  conclusion  that  1/2  of  all 
collisions  of  Nf^D,  ^P]  +  O2  lead  to  NO  with  little  or  no  rovibrational  excitation.  Detailed 
studies  of  these  reaction  surfaces  are  required  to  further  analyze  the  branching  rates  in  this 
system. 
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SECTION  5 


RECOMMENDATIONS 


Previous  studies  of  the  electronic  structure  and  energy  levels  of  the  uranium/oxygen  system 
have  been  undertaken  by  UTRC  and  reported  in  DNA-TR-82-159  (Reference  85), 
DNA-TR-85-156  (Reference  86),  DNA-TR-88-12  (Reference  87)  and  the  present  technical  report 
for  Contract  DNA(X)l-88-€032.  In  these  studies  the  focus  was  on  estimating  the  LWIR  radiation 
characteristics  and  the  efficiency  of  solar  pumping  of  the  allowed  electronic  transitions  in  this 
system  which  fall  in  the  visible  wavelength  region.  An  examination  of  the  thermodynamics  and 
kinetics  of  the  uranium/oxygen  system  shows  that  U'*'  and  UO'*'  are  terminal  ions,  and  therefore 
potentially  important  radiators  at  high  altitudes,  whereas  UO2'''  is  more  important  at  low 
altitudes  where  there  is  significant  O2  concentration.  Further  studies  of  the  uranium/oxygen 
system  are  outlined  below.  In  addition,  suggested  studies  of  dielectronic  recombination  and 
kinetic  branching  in  atmospheric  reactions  producing  radiating  NO  molecules  are  described. 

5.1  URANIUM/OXYGEN. 

The  calculations  of  the  electronic  structure  of  the  uranium  oxide  ions  have  been  performed 
with  neglect  of  the  spin-orbit  interaction  term,  optimizing  the  relativistic  density  functional  SCF 
wavefunctions  for  the  particular  ion.  Spin-orbit  interaction  is  then  introduced  as  a  perturbation 
through  configurational  mixing.  This  same  approach  has  been  used  by  Krauss  and  Stevens 
(Reference  89),  with  a  full  initio  treatment  of  the  valence  electrons,  using  /-dependent 

relativistic  effective  core  potentials.  Either  approach  allows  for  polarization  and  charge  transfer 
from  the  ligand  oxygen  atoms  to  the  central  uranium  ion. 

A  different  approach  for  rare  earth  metal  coddes  has  been  described  by  Dulick  (Reference 
105)  which  is  based  on  a  crystal  field  model.  Application  of  this  method  to  the  UO  and  UO'^ 
ground  states  (Reference  106)  yields  a  significantly  different  energy  ordering  for  the  Q  coupled 
states  than  is  found  in  the  full  molecular  treatments  of  the  systems  described  above.  The 
application  of  the  crystal  field  model  to  lanthanide  monoxides,  however,  appears  to  yield  the 
known  experimental  electronic  assignments  (Reference  107).  An  examination  of  the  composition  of 
the  predicted  Q  states  for  UO'*',  based  on  this  crystal  field  model,  indicates  a  significant  mixing 
of  electronic  states  of  the  uranium  ion  of  small  /-value.  These  efforts  are  not  found  in  the 
molecular  treatments,  whose  compositions  are  dominated  by  the  tightly  coupled  5f  states.  Further 
studies  of  these  differences  in  approach  are  indicated;  however,  since  all  methods  are  either  based 
on  or  predict  a  large  ionic  character  in  the  UO'*'  and  UO2'''  ions,  our  conclusion  of  strong  LWIR 
emission  should  not  be  altered.  The  crystal  field  model  has  not  been  utilized  for  describing 
electronic  transitions  or  the  absorption/emission  characteristics  in  either  the  IR  or  visible 
wavelength  regions. 
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Our  suggested  studies  include  the  further  examination  of  low-lying  electronic  states  in  UO 
and  UO2''’.  These  studies  will  primarily  cover  the  IR,  visible  and  UV  excitation  wavelengths 
( <  -  3  eV)  in  order  to  more  accurately  define  the  radiation  emission  strengths.  A  new  density 
functional  computer  code,  has  been  developed  in  collaboration  with  Professor  F.  E.  Harris 
(University  of  Utah),  and  would  be  utilized  to  carry  out  these  suggested  studies.  This  new  code 
yields  a  more  accurate  representation  of  the  electron  density  in  the  valence  region  between  atomic 
centers.  This  region  has  previously  been  described  by  a  constant  (muffin-tin)  potential  which  can 
lead  to  inaccuracies  in  describing  polarization  effects.  We  feel  that  this  new  code  would  also 
permit  a  more  critical  evaluation  of  the  crystal  field  model,  as  applied  to  the  uranium/oxygen 
system. 

Our  increase  in  the  data  base  for  the  uranium/oxygen  system  will  include  radiation 
f-numbers,  characteristic  absorption/emission  wavelengths  and  calculated  density  of  states  for 
U"*"",  UO"^",  U02''‘"  (n  =  0,1,2).  These  studies  would  be  aided  by  a  parallel  experimental 
program  on  the  spectroscopy  of  electronic  states  of  UO  and/or  UO'*'.  A  suggested  collaborator  is 
Professor  Michael  Heaven  at  Emory  University. 

5,2  DIELECTRONIC  RECOMBINATION. 

In  addition  to  our  suggested  calculations  of  the  structure  and  radiative  properties  of 
uranium/oxygen  species,  further  analysis  of  the  importance  of  dielectronic  recombination 
processes  should  be  carried  out.  A  preliminary  study  of  the  important  dissociative  recombination 
states  of  U°,  U'*'  and  U'*"*'  has  been  described  above  under  Section  4.  Several  highly  excited 
configurations  for  U°  have  been  identified  which  indicate  that  e  +  U'*'  may  be  an  important 
charge  neutralization  process,  but  a  more  extensive  study  is  required  to  determine  the  density  of 
states  accessible  to  dielectronic  recombination  under  conditions  up  to  perhaps  SOGO^K. 

A  corresponding  DR  process  is  possible  for  e  +  UO'*'.  At  the  present  time,  only  the 
optically  connected  excited  states  of  UO  have  been  studied.  The  ground  molecular  orbital 
configuration  for  UO"^  is: 


l>n  \<p  \6]  (114) 

where  the  last  three  molecular  orbits  are  nearly  pure  f-orbitals  centered  on  uranium.  Excitation 
processes  that  need  to  be  examined,  all  in  a  molecular  framework,  are:  -*  d^, 

^  -•  s^,  f^  -►  sd,  f^  -•  p^,  f^  -•  sp,  and  f^  -*  pd.  Additional  excitations  from  the  core  2jt  level 
will  be  examined  but  it  is  expected  that  these  levels  are  too  tightly  bound  for  the  electron  energy 
range  of  interest. 
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These  studies  of  dielectronic  recombination  contribute  to  our  understanding  of  the  kinetics 
of  charge  neutralization  following  nuclear  bursts.  In  particular,  the  role  of  dielectronic 
recombination  as  a  mechanism  for  removal  of  electrons  is  uncertain  at  present  and,  if  found  to  be 
important,  must  be  included  in  the  chemistry  kinetics  in  the  NORSE  code. 

S3  NO/NO  +  FORMATION. 

Under  this  Contract,  DNA001-8S-C-0032,  a  study  of  the  atmospheric  reactions  O  4-  N2  and 
N  +  O2  which  yield  vibrationally  excited  NO  molecules  was  initiated.  This  study  was  in  response 
to  the  requirement  defined  by  Gilmore  (Reference  108)  for  a  more  accurate  analysis  of  the 
NpD]  +  O2  and  NpP]  +  O2  reactions  which  are  modeled  in  the  NORSE  calculation  of  NO 
chemiluminescence.  Since  the  N[^P]  +  O2  reaction  can  yield  NpD]  +  O2  as  products,  as  well  as 
vibrationally  excited  NO,  there  is  a  large  uncertainty  in  the  kinetic  branching.  If  vibrationally  hot 
NO  is  strongly  favored  as  the  product  channel,  high  IR  emission  at  9/<m  and  4.5  ;/m  is  predicted. 
Our  initial  studies  of  the  kinetics  of  the  N[^]  +  O2  reaction  suggest  that  only  1/6  of  all  collisions 
involving  N[^P]  atoms  appear  to  lead  to  rovibrationally  excited  NO  molecules.  Thus  the  present 
atmospheric  codes  predict  too  much  excitation  for  the  NO  product  of  this  reaction. 

In  a  closely  related  problem,  a  strong  radiator  in  the  IR  following  a  nuclear  burst  is 
vibrationally  excited  NO'*',  with  emission  near  43 //m.  An  examination  of  the  ion-molecule 
reaction  surface  for  O  +  N2'*’,  Oi"^  +  N,  and  O'*’  +  N2  should  be  carried  out  to  establish  the 
overall  product  correlations  for  this  system.  If  required,  detailed  quantum  mechanical  studies  of 
the  pertinent  reaction  surfaces  should  be  examined.  A  similar  study  of  the  O'*"  +  N2-*  NO^  + 

N  reaction  has  been  previously  carried  out  for  DNA  by  UTRC  under  AFGL  Contract 
F19628-80-C-0209  (Reference  109).  All  of  these  studies  of  atmospheric  molecules  have  the 
objective  of  increasing  the  fundamental  data  base  of  the  reaction  kinetics  in  order  to  provide 
useful  input  for  chemiluminescence  modeling  in  the  various  atmospheric  codes. 
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Figxire  1.  Potential  energy  curves  fw  UO+, 
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Allen  etal.:  JCP.  89.  p.  5363,  1988 


Figure  2.  He  I  photoclectron  spectrum  of  the  vapor  above  an  equimol 
mixture  heated  to  2250  ±  30  K. 


a)  b)  c) 


a)  Incoming  electron  has  energy  less  than 

b)  Incoming  electron  4-  KE  in  Coulomb  field  excites  an  electron 
from  nl  — >  nl'  and  incident  electron  is  trapped  in  a  highly 
excited  state. 

c)  Doubly  excited  state  radiates  (nl'  >  nl)  resulting  in  a 
stable  excited  state  lying  below  the  ionization  threshold. 


Figure  3.  Dieleccnonic  recombination. 
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Figure  4.  Molecular  correlation  diagram  for  low-lying  states  of  NO2  in  symmetry. 
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POTENTIAL  ENERGY.  E  (eV) 


Potential  energy,  eV 


Figure  6.  Molecular  correlatitm  diagram  for  low-lying  doublet  states  of  NO2  in  Cj  symmetry. 
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Table  1.  Calculated  oscillator  strengths  for  the  -  l^K  transition  of  UO^. 
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Table  2.  Calculated  oscillator  strengths  for  the  Hi  -  3‘*tl  transition  of  1)0'''. 
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Table  3.  Calculated  oscillator  strengths  ^v  v*)  fo**  Ihc  I**!  -  2‘*H  transition  of  UO'*'. 
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Table  5.  Calculated  oscillator  strengths  f^v  v")  -  2*1  transition  of  UO'*'. 
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Table  6.  Calculated  oscillator  strengths  r(v<v»)  for  the  1^1  -  1^11  transition  of  CO^. 
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Table  7.  Calculated  oscillator  strengths  f^v-v")  for  the  1^1  -  3^1  transition  of  DO'*'. 
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Table  8.  Calculated  oscillator  strengths  l(v  v»)  for  the  l^H  - 1^\  transition  of  UO^ . 
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Table  9.  Calculated  oscillator  strengths  Ibc  -  2^H  transition  of  UO'*'. 
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Table  10.  Calculaled  oscillator  strengths  f(v'v')  for  the  I'^H  -  l^K  transition  of  UO'*'. 
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Table  11.  Calculated  oscillator  strengths  f(v>v'')  for  llx  l^H  -  l^H  transition  of  DO'*'. 
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Table  13.  Calculated  oscillator  strengths  f(v>v")  the  -  3^1  transition  of  liO^. 
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oscillator  strength 
wavelength,  microns 


Table  14.  Energetics  of  dielectronic  recombination. 


e  +  U+ 

€  (cm”*) 

U* 

=  +  [5fi/a7s2] 

0. 

[5f^/2'7s7pi/27p3/2] 

-0. 

(5ff/27s^7p3/26d5/2] 

116. 

l5fi/27s7pi/2] 

6154. 

[5f|/27s6d3/2n/l 

9589. -E„, 

[5f|/27s6d5/2n/] 

14071.  - 

[5ff/27s^3/2nfl 

• 

19924. -E„, 

e  +  U+  + 

€  (cm'^) 

[u^r 

€  (cm'^) 

e  +  [5fi/a7sl 

0. 

{5^^27s7p3/2] 

2711. 

I5^3/26d3/2nfI 

4117. -E„, 

[5f|/26ds/2nf) 

9345. -E,rf 

[5^/27s6d3/2nf) 

24660. -E,^ 

[5^/2^3/2"fl 

2472Z-Eb, 

!  +  I5f5/26d3/2] 

4117. 

[5fi/26d5/2nfl 

5228. -E,rf 

[5fi/27s6d3/2nfl 

20543. -E,rf 

[5f|/26di/2n^ 

20604.-E,rf 
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Table  IS.  e  +  U'*'  dielectronic  recombination  rate  (preliminary  data). 


ot  (cm^/sec) 


T(K) 

5f37s7p2 

5f27s27p6d 

300 

5.9  X  10-11 

1.6  X  10-23 

1000 

1.4  X  10-11 

2.4  X  10-13 

10000 

5.1  X  10-13 

2.2  X  10-13 

20000 

1.8  X  10-13 

3.4  X  10-13 
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Table  16.  Molecular  correlation  diagram  for  N  +  O2. 


Reactants 

N('Su)  +  O2 

N('^)  +  02(a^Ag) 

N(%)  +  02(bi2g*) 
N(2Du)  +  02(X32^ 

N(2Du)  +  02(a^Ag) 
N(2Pu)  +  02(X32^ 

N(2Du)  +  02(bi2g+) 

N(-»Su)  +  02(c^tj 

N(^,)  +  02(A'  3Au) 

N(%)  +  02(A32;u+) 


Energy  (eV) 

States 

Point  Group 

1.401 

2.4.62+ 

Coov 

2.4.6B1 

Czv 

2,4.6a< 

Cs 

1383 

Cqov 

Czv 

^A",  '*A' 

Cs 

3.037 

42- 

Coov 

^Az 

Czv 

4a» 

Cs 

3.785 

2,42  +  ,  2.4n,  2.4^ 

Coov 

2.^Bi  (2),  2.4B2,  2.4A2.  2.4Ai, 

Czv 

2.^A"  (2),  2.4a'  (3) 

Cs 

4.766 

22+,  22-,  2n.  2a,  2<I),  2r 

Coov 

2Ai(2),2A2(3).2Bi(3),2B2(2X 

Czv 

2A'  (5),  2a"  (5) 

Cs 

4.977 

2,42-,  2,4n 

Coov 

2.4Ai.2.4A2,2.4B2 

Czv 

2-4A',  2.4A''  (2) 

Cs 

5.420 

22-,  2a  2a 

Coov 

2Ai,2A2(2),2Bi,2B2 

Czv 

2A'  (2X  2A"  (3) 

Cs 

5.499 

42  + 

Coov 

4Ai 

Czv 

4A' 

Cs 

5.701 

2,4.6a 

Coov 

2.4.6A1,  2.4.6B2 

Cjv 

Z4.6a'  2,4,6a  " 

Cs 

5.789 

2.4.62- 

Coov 

2.4.6B2 

Czv 

i4,6A» 

C, 
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Table  16.  Molecular  correlation  diagram  for  N  +  O2.  (Concluded) 


Reactants 

Energy  (eV) 

States 

Point  Group 

N(^Pu)  +  02(aiAg) 

5.958 

2n.  2a.  2<|) 

Coov 

2Ai(2),2A2.2Bi,2B2(2) 

C2v 

2A'  (3).  2A''  (3) 

Cs 

N(2Pu)  +  02(bi2g+) 

6.612 

2s+.2n 

Coov 

2Ai.  2Bi,  2B2 

C2V 

2A'  (2),  2A" 

Cs 

Products 

Energy  (eV) 

States 

Point  Group 

NO(^  +  0(^Pg) 

0 

2,42+  2,42-  2,411^  2,4^ 

n 

^oov 

2>^A'(3X  ^^A''(3) 

Cs 

NO(2n)  +  0  (iDg) 

1.967 

22  +  ,22-  2n(2X2A,2<I> 

Coov 

2A'  (5X  2a-  (5) 

Cs 

NCK^TI)  +  0  OSg) 

4.189 

2n 

Coov 

2A'.  2a- 

Cs 

NO  (a^n)  +  0  (^Pg) 

4.765 

2,4,62  +  ^  2,4,62-  2,4,6fl^  2,4,6^ 

Coov 

2’'^’^A'  (3X  2.4.6^'  (3) 

Cs 

N0(A22+)  +  0(3Pg) 

5.450 

2,42-.  2,4n 

Coov 

2-^A'.  2.^A-  (2) 

Cs 

NO  (B^n)  +  0  (3Pg) 

5.693 

2,42  +  ^  2,42-  2,4n,  2,4^ 

Coov 

2-^A'  (3),  2.*A-  (3) 

Cs 

NO  (b^2-)  +  0  (3Pg) 

6.0344 

2,4,62  +  ,  2,4,6n 

Coov 

2'‘»’6A'  (2),  2-4.6a» 

Cs 
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Table  17.  Adiabatic  reaction  surfaces  arising  frtkm  N(^P]  +  O2. 


N  *  [2p]  +  02[X323 


N02*(^A'X]  _ 

NO2  •  f^A'  VI)  - 
NO2  •  [^A"  vin)  -► 

N02*[^A"V]  ^ 

N02*pA"IX]  - 
N02*[^A'V1]  - 


NO  •  [a.^n\  +  O  [^P]  +  0.212  eV 
NO  •  [a^ri]  +  O  I^P]  +  0.212  eV 
NO  [X^ri]  +  0*(‘D]+  3.010  eV 
NO  •  [a‘*n]  +  O  [3p]  +  0.212  eV 
NO  [X^n]  +  0*[^S)  +  0.788  eV 
NO*(a^n]  +0  [3P]  +  0.212  eV 


Tabie  18.  Adiabatic  reaction  surfaces  arising  firom  N[^D]  +  02> 


N«[^l  +  02(X32j 


n02*[^A'  n) 
N02*pA'  ffl] 
N02*[2A'  ivj 

N02*pA'  I] 

n02*[^A”  n] 

N02*[^A'  ID] 
N02*[^A'  IV] 
N02*[^A'  V] 

n02*[^A"  ni] 

NOi’l^A'  IV] 


NO  [x^n]  +  o 
NO  pc^ri]  +  o 
NO  (x^n]  +  o 

NO  pC^Tl]  +  O 
NO  PC^TI]  +  O 
NO  PC^TI]  +  O 
NO  •  [a^n]  +  O 
NO*[a^n]  +  O 
NO  [X^Tl]  +  O 
NO*[a^n]  +  O 


[3P]  +  3.785  eV 
[3p]  +  3.785  eV 
{1D]+  1.878  eV 
[3p]  +  3.785  eV 
[3p]  +  3.785  eV 
[3p]  +  3.785  eV 
[3p]  -  0.980  eV 
[3P]  -  0.980  eV 
[3p]  +  3.785  eV 
[3p]  -  0.980  eV 
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Table  19.  Minimum  energy  long  range  interaction  potentiais  for 
N  +  O2  in  Cj  symmetryi 


2A' 

R  =  4.0  bohrs* 

R  =  3.0  bohrs 

Long  Range 
Character 

%  +  X^Sg 

-201.4657103 

-201.4657669 

a 

-201.3423547 

-201.3428036 

a 

2Du  +  X^Zg 

-201.3423540 

-201.342542 

r 

-201.3423535 

-2013415228 

r 

-201.3045858 

-201.3051402 

a 

-201.3045827 

-2013048667 

a 

2Du  +  a^Ag 

-2013045808 

-2013043220 

r 

-2013045790 

-2013041826 

r 

-201.3045775 

-201.3040453 

r 

2Pu  +  X^Zg 

-201.2994959 

-201.2988189 

r 

^A" 

2Da  +  X^Zg 

-201.3423546 

-201.3427958 

a 

-2013423536 

-201.3420978 

r 

-2013045845 

-2013051136 

a 

-2013045830 

-2013048777 

a 

^Du  +  a^Ag 

-2013045810 

-201.3043588 

r 

-201.3045790 

-2013041614 

r 

-201.3045775 

-2013040561 

a 

2Pu  +  x^z; 

-201.2994986 

-201.2999683 

a 

-201.2994921 

-201.2987603 

r 

®R  is  measured  from  N  to  the  center  of  the  O2  molecule,  energy  in  hartrees. 
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